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Abstract 


This article reviews our understanding of the characteristics and causes of northern Eurasian 
summertime heat waves and droughts. Additional insights into the nature of temperature and 
precipitation variability in Eurasia on monthly to decadal time scales and into the causes and 
predictability of the most extreme events are gained from the latest generation of reanalyses and from 
supplemental simulations with the NASA GEOS-5 AGCM. Key new results are: 1) the identification 
of the important role of summertime stationary Rossby waves in the development of the leading 
patterns of monthly Eurasian surface temperature and precipitation variability (including the 
development of extreme events such as the 2010 Russian heat wave), 2) an assessment of the mean 
temperature and precipitation changes that have occurred over northern Eurasia in the last three 
decades and their connections to decadal variability and global trends in SST, and 3) the quantification 
(via a case study) of the predictability of the most extreme simulated heat wave/drought events, with 
some focus on the role of soil moisture in the development and maintenance of such events. A 
literature survey indicates a general consensus that the future holds an enhanced probability of heat 
waves across northern Eurasia, while there is less agreement regarding future drought, reflecting a 
greater uncertainty in soil moisture and precipitation projections. Substantial uncertainties remain in 
our understanding of heat waves and drought, including the nature of the interactions between the 
short-term atmospheric variability associated with such extremes and the longer-term variability and 
trends associated with soil moisture feedbacks, SST anomalies, and an overall warming world. 
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1 . Introduction 


“ While in western Europe there is continual rain and they complain about the cold summer, here in 
Russia there is a terrible drought. In southern Russia all the cereal and fruit crops have died, and 
around St Petersburg the forest fires are such that in the city itself, especially in the evening, there is a 
thick haze of smoke and a smell of burning. Yesterday, the burning woods and peat bogs threatened 
the ammunition stores of the artillery range and even the Okhtensk gunpowder factory P 2 . This 
remarkable July 15, 1875 entry in General Dmitry Milyutin’s diary reflects not only the fact that 
Russia suffered from terrible drought and heat in the past, but also a realization long ago that such 
droughts were at times juxtaposed with cool and wet conditions over Europe. Today, we know this 
juxtaposition is no coincidence but in fact reflects the unique large-scale atmospheric controls on 
drought and heat waves affecting much of northern Eurasia. Droughts in Eurasia indeed have a 
character all their own. 

Historical records show that over time the peoples of the Eurasian continent have suffered through 
numerous heat waves and droughts, events that have impacted the course of battles 3 , desiccated 
important crop lands (thereby inducing famine), produced numerous forest and peat fires, and 
contributed to thousands of deaths. Gumilev (1960) and Pines (2012) review the pulses of dry and 
relatively humid periods that have occurred during the past 2 millennia over the entire Great Steppe of 
northern Eurasia (from Pannonia in the west to Manchzhuria in the east), causing prosperity and decay 
of ancient states and the migration of nomadic tribes. A compendium of extreme events during the 
past 1100 years over ancient European Russia (ER), Belarus, and Ukraine (so-called Kiev Rus’) is 

2 http ://there se-phil .li vej oumal .com/171196 .html 

3 The town of Szigetvar, Hungary was under siege by the Turks on 7 August 1566. The main protection of the town was a lake 
and marshland that normally surrounded it. “Chance however now favored the Turks. A drought had prevailed during the two 
preceding months, and the terrain surrounding the old town had become so dry, as considerably to facilitate the approach of the 
enemy.” Vambery (1886). 
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provided by Borisenkov and Pasetsky (1988), Bogolepov (1922), and Vazhov (1961). Appendix B 
contains a list of the major droughts and heat waves that have occurred since the late 19 th century. 

Droughts continue to have major impacts on northern Eurasian agriculture. As noted in Golubev and 
Dronin (2004), “Another notable feature of Russian agriculture are the rather large fluctuations in 
year-to-year yield, which are considerably higher than in any other major grain producing country in 
the world... These high fluctuations in total cereal production were undoubtedly the result of irregular 
precipitation.” In fact, many of the important early studies on Russian drought and temperature 
extremes were performed to address their impacts on agriculture in the important growing regions of 
Povolzhie, North-Caucasus and Central-Chernozem, regions that produce about 2/3 of the Russian 
food grains (Kleshenko et al. 2005). Kahan (1989) lists some major Russian droughts (based on the 
work of Rudenko 1958; see also Appendix B) and their agricultural impacts; he notes that the 
increased impact of natural calamities was in part associated with the expansion of Russian agriculture 
(mainly grain acreage) toward the south and southeast into the steppes and semi-arid regions 
characterized by drier climatic conditions. However, these regions are also characterized by fertile 
soils (chernozem) and have longer growing seasons 4 . This, together with a general increase of grain 
productivity, appears to have made the severe drought-induced famines of previous years much less 
likely. For example, the 2010 heat wave and drought over ER (which was so severe that we had to go 
back to 1092 AD to find an analogue) 5 caused not a famine but did cause a stoppage of grain export 


4 For example, when at the end of the rule of Tsar Boris Godunov the Moscow Tsarstvo was struck by a 3-year long famine (1601- 
1603), cold summers were to blame, most probably related to a catastrophic volcanic eruption of Huaynaputina (Peru) in 1600. At 
that time the present major grain areas of ER were not plowed and the present region of “sustainable agriculture” in the forested 
areas of central ER around 55°N simply did not receive summer temperatures sufficient for grain harvests. Ultimately, this famine 
caused an 8-year-long period of social turmoil, civil war, and invasion by Swedish and Polish marauders, and it ultimately caused 
a change in the ruling dynasty. 

5 National Yearbooks (Letopisi) prepared by Russian monks since the early 10 th century report most important political and 
environmental events over Kiev Rus’ (the area of present northern Ukraine, eastern Belarus, and central part of European Russia). 
In the 20 th century, Letopisi were summed and reanalyzed (Borisenkov and Pasetsky 1988; Barash 1989). Generally, summers 
over Kiev Rus’ in the 1 1 th century were mostly warm and dry. However, on this background, the 1092 summer was extremely dry. 


4 


77 

78 

79 

80 

81 

82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 


from Russia. Multiple-year droughts have occurred over the past century, especially in the heartland 
of the Eurasian steppes (Kahan 1989). Such droughts include the 5-year event (1929-1933) in the 
Akmolinsk (presently Astana, the capital of Kazakhstan) area. Historically, famines, when they occur, 
tend to be associated with multiple years of drought (Borisenkov and Pasetsky 1988). 

An interesting and telling aspect of the literature addressing droughts in northern Eurasia is the lack 
therein of a clear distinction between drought and heat waves. To some extent, this is because 
summer dryness in this region has two different manifestations: agricultural drought (i.e., soil 
moisture deficits) and “fire weather” (in the forested areas of northern Eurasia, a prolonged period of 
hot weather with little or no rainfall) as described in, for example, Nesterov (1949) and Groisman et 
al. (2007). Agricultural droughts in northern Eurasia also may last for several weeks or even months, 
particularly under conditions of a short growing season where a complete harvest loss can be caused 
by a short heat wave that strikes at a critical period of wheat development. Another aspect of droughts 
in the steppe and semi-desert zones of northern Eurasia are “sukhovey” - an extended period of dry 
hot winds characterized by intense transpiration and rapid wilting of vegetation (Lydolph 1964). 
Sukhovey typically emanate from the periphery of anticyclones, bringing in warm and dry air 
originating in the deserts of Africa, Asia Minor and southern Kazakhstan 
(http://en.wikipedia.org/wiki/Sukhovey) . Historically, sukhovey have been a major impediment to 
large-scale sedentary agriculture in Central Asia (Sinor 1994). An important point here is that the 
traditional notions of meteorological (precipitation deficits) and agricultural droughts (soil moisture 

Moscow Letopis summary for 1092 says: “Huge circle was in the sky in this summer, a drought was so strong that soil was burned 
and many forest and swamps were set in fire themselves”. Letopisi witness: clear skies throughout the entire summer; prolonged 
period without rainfall; extremely hot weather; fields and pasture “fired out “, and widespread naturally caused forest and peat bog 
fires (let us recall that at that time wetlands were undisturbed which is opposite to the present state of affairs). In Kiev, in the 
following autumn and winter more than 7 thousand (of total 50,000) died from starvation. Losses beyond the capital city were (in 
percent) even higher. This unfortunate development was followed by widespread epidemics. 
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deficits) are perhaps not as relevant nor as clearly separated in northern Eurasia as in other regions of 
the world. 


The strong link between heat waves and drought in northern Eurasia suggests that we should treat 
them as different facets of the same phenomena. In fact, many metrics of drought in northern Eurasia 
involve an explicit temperature criterion; a drought is said to occur, for example, only after a certain 
minimum number of days with temperatures above a certain threshold (e.g., Selianinov 1928). The 
strong connection between drought and excessive heat reflects in part the central role of anticyclones 
in the development of northern Eurasian droughts; the anticyclone inhibits precipitation by blocking 
or diverting the westerlies and storm systems, and it increases temperature through descending 
motions (which further inhibit precipitation) and increased insolation associated with clear skies (e.g., 
Buchinsky 1976). Another relevant mechanism involves soil moisture feedback on temperature; 
reduced precipitation leads to reduced evaporative cooling of the land surface. While multi-year 
droughts do occur in Eurasia, particularly toward the south of our study area, most droughts have 
shorter time scales; most severe events 6 occurring across northern Eurasia in fact rarely exceed 50 
days in duration (Cherenkova, 2007). The north/south differences and east/west differences in 
drought occurrence reflect the spatially varying influences of the oceans and various air masses 
(tropical, subtropical, polar) across the continent and the arrangement of Eurasia’s major mountain 
chains. 

The recent extreme heat waves and droughts of 2003 (Europe) and 2010 (Russia) have highlighted the 
urgency of understanding better their causes and whether or not they are a manifestation of a wanning 

6 However, as mentioned above , it is enough in this part of the world to have several days of adverse weather during particular 
periods of cereal growth to cause crop loss. Such short-term events (waves) are much more frequent than long periods of 
decremented weather. 
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world (e.g., Dole et al. 2011; Trenberth and Fasullo 2012; Otto et al. 2012). While there is 
agreement on many aspects of such extreme events, including the role of anticyclones, there are still 
substantial unknowns about their causes and predictability. In particular, we do not yet understand 
which large scale processes (including climate change, SST, monsoons, links to higher latitudes) may 
have played a role in making them so exceptional. 

In this paper, we delve into some of the outstanding questions regarding the nature and causes of 
Eurasian heat waves and droughts, their predictability, and what we can expect in a future warmer 
world. We focus on northern Eurasia, in particular the region outlined in the Northern Eurasia Earth 
Science Partnership Initiative (NEESPI, Groisman et al. 2009) - longitudinally, this region extends 
from 15°E to the Pacific coast, and latitudinally it extends from 40N to the Arctic Ocean coastal zone. 
The region includes the territory of the fonner Soviet Union, Fennoscandia, Eastern Europe, Mongolia 
and north China, though our presentation of the results often extends beyond it in order to provide a 
more global perspective of relevant teleconnections and physical mechanisms. 

This paper is part of a Global Drought Information System (GDIS) special collection that addresses 
the causes of drought world-wide. We note that there are separate papers in this collection focusing 
on drought in large regions bordering and in part overlapping northern Eurasia, including papers on 
Europe, the Middle East and southwest Asia, and eastern Asia (Seneveratne et al. 2013; Barlow et al. 
2013; Zhang and Zhou 2013). The interested reader is referred to those papers for more infonnation 
on those regions. In keeping with the guidelines of the submissions to the GDIS special collection, 
we will touch on the following topics: 

i) Drought/heat wave occurrence, metrics, and impacts 
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ii) Key regional circulation features and physical processes 

iii) Trends 

iv) Predictability and projections 

v) Gaps in our understanding 

We begin in Section 2 by discussing the morphology and metrics of droughts and heat waves. Section 
3 examines the physical mechanisms responsible for their occurrence. Section 4 contains a review 
and analysis of interannual variability and trends, and Section 5 discusses projections and 
predictability. A summary is provided in section 6. In addition, two appendices are provided: 
Appendix A describes the datasets and model simulations used in this study, and Appendix B provides 
a compilation (based on various sources including research papers and the popular literature) of some 
of the major droughts and heat waves that have occurred in northern Eurasia since 1875. 

2. The morphology and metrics of northern Eurasian heat waves and droughts 

a) Characterizing drought and heat waves 

The previous discussion highlighted the adverse impacts of prolonged drought conditions, as well as 
shorter-period (weeks to months) heat waves and related droughts, on the main agriculture regions of 
northern Eurasia, with the latter events also playing an important role in the occurrence of “fire 
weather” in the forested areas of northern Eurasia. The discussion also emphasized the importance of 
the interplay between temperature and precipitation variability in the development of droughts. 
Droughts are ultimately driven by precipitation deficits, and their impacts (e.g., on agriculture) depend 
on the extent to which they lead to deficits in soil moisture and other water resources important to 
society. Temperature increases associated with the precipitation reductions, which can act to 
exacerbate the drought conditions, can result from reductions in cloudiness; reduced cloud cover over 
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northern Eurasia in summer can warm the surface, since daytime warming is 2-3 times stronger than 
the nighttime cooling, and the period of daylight is longer than the period of nighttime (Groisman et 
al. 1999; Tang et al. 2012; Tang and Leng 2012a, b). The temperature increases can also result from 
reductions in soil moisture, which lead to reduced evaporation and thus reduced evaporative cooling 
of the surface. Finally, various dynamical processes linked to (for example) the development of 
persistent anticyclones (the subject of the next section) can act to further reduce cloudiness and 
precipitation and at times can lead to intense heat waves. 

Here we look at various measures that highlight and quantify the above aspects of droughts and heat 
waves over northern Eurasia. A number of different indices are available for consideration. Most are 
based on joint consideration of precipitation and temperature, with a focus on agricultural applications 
(e.g., Meshcherskaya and Blazhevich 1997; Kleshenko et al. 2005); examples include the 
hydrothermal coefficient (HTC, Selianinov 1928) and the dryness index of Ped (1975). 

Meshcherskaya and Blazhevich (1997) developed a combined drought and excessive moisture index 
(DM) that takes into account the areal extent of the precipitation and temperature anomalies. 

Appendix B provides more details on these and several other of the more popular metrics. 

The overall character of northern Eurasian precipitation and temperature variability and its 
relationship to drought and heat waves is now examined. Figure 1 (left panels) shows the variance of 
June -July-August (JJA) mean precipitation for the last three decades (1979-2012) as determined by 
two different reanalyses and from an analysis of station observations. While there are some 
differences, there is general agreement that the largest precipitation variance over northern Eurasia 
occurs over European Russia extending eastward along about 55°N - a region for which the maximum 
rainfall in summer is associated with cyclones from the Atlantic reaching the Yenissey River Valley in 
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Central Siberia. Other regions with relatively large precipitation variance are found in the West 
Caucasus and in the mountainous regions of southeastern Russia/northeast China. Minimum 
variances are seen in the desert regions east of the Caspian Sea extending to northwestern China and 
Mongolia. 

While a precipitation deficit is the basic ingredient of drought, the link between the magnitude of such 
deficits and the presence of drought, at least as defined by established drought indices, is not readily 
apparent. To make this link quantitative we take as an example the connection of the precipitation 
deficits with the Ped (1975) dryness index, S, (the difference between the nonnalized anomalies of 
JJA-mean surface temperature and precipitation - see Appendix B for definitions). The variance of 
Si is 

V = 2 (l- p) , (1) 

where p is the correlation between the surface temperature and precipitation. The variance in drought 
as measured by the Ped (1975) dryness index thus highlights the fact that drought depends on the 
interplay between precipitation and temperature. Figure 1 compares the variance field for the drought 
index (right panels) to that for precipitation (left panels). Note that a comparison of the magnitudes is 
irrelevant, since the fields have different units; here, we consider only the comparison of spatial 
patterns. For the dryness index, the largest values (regions with the largest negative correlations 
between P and T) occur farther south, extending across the main agricultural regions of northern 
Eurasia, as well as over much of Mongolia and northeast China. Consistent across the three estimates 
are the relatively large values to the north and east of the Caspian Sea (including northern 
Kazakhstan), over the Caucasus, and over much of Europe, including the Balkans. In essence, the 
comparison in Figure 1 shows that standard drought indices do not simply describe the absolute 
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magnitudes of precipitation deficits but, as we shall see in Fig. 2, also reflect the important link with 


temperature. 

We note that the coefficient of variation (C v ) relates precipitation fluctuations to mean or “normal” 
climate values - a metric that is in principle 7 perhaps a better measure of precipitation irregularity as it 
relates to potential impacts on agriculture. Shver (1976) showed, for example, that while C v of 
monthly rainfall totals over Eurasian agricultural areas north of 55 "N from May to July are close to 
0.5, further south C v gradually increases and in the North Caucasus, southern Ukraine, and northern 
Kazakhstan reaches values of 0.7 in July and 0.8 in August, indicating the latter regions experience 
larger swings in precipitation relative to the mean. 

We next turn to the temperature variance. The left panels of Figure 2 show that the variance of the 
JJA mean is characterized by generally increasing values with latitude, with the largest variance 
occurring north of Caspian Sea and over the Ural Mountains. While this looks nothing like the 
distribution of the precipitation variance (Fig. 1, left panels), we have already seen that there are 
regions where the precipitation and temperature are correlated (Fig. 1, right panels). We can quantify 
the extent to which the temperature variability is “explained” by precipitation variability via simple 
linear regression. The results (left panels of Figure 2) indicate that precipitation variability explains a 
substantial fraction of the temperature variance over much of southern European Russia and western 
Siberia (e.g. through evaporative cooling). The similarity to the variance of the Ped index (c.f. Fig. 1, 
right panels) is not surprising since both measures depend on the correlation between the temperature 
and precipitation. Furthermore these regions occur in the transition between so-called water-limited 

7 In practice this metric is sensitive to bias in the estimates of the mean state. 
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(to the south) and energy limited (to the north) climate regimes (Koster et al. 2006a, Fig. 4a), where 
land feedbacks are particularly important (Koster et al. 2004). We will look more directly at the 
impact of soil moisture feedbacks on temperature variability in the context of model simulations in 
section 3. 

Of course, as mentioned above, cloud cover can also affect temperature variability. In fact, in an 
assessment of the role of cloud cover and rainfall on the daytime temperature (T max ), Tang et al. 
(2012) showed that for the western half of northern Eurasia (where the major agricultural regions 
reside), summer cloud cover is negatively correlated with T max and that these correlations are much 
stronger than those with precipitation. Tang and Leng (2012b) show that the variance of Eurasian 
summer T max is better explained by changes in cloud cover than by changes in precipitation at high 
latitudes and in the middle latitude semi-humid area, while in northern Eurasia the dependence on 
precipitation is strong only in the Central Asia arid area. 

The above findings suggest that heat waves in Northern Eurasia are influenced by both soil moisture 
(and precipitation) and circulation (and cloud cover) anomalies, though it is still unclear which plays 
the more important role. The interactions involved are indeed complex; precipitation deficits can be 
caused by decreases in cloudiness, and a dry land surface can suppress evapotranspiration and thus 
inhibit local cloud formation. All said, it seems reasonable to pay significant attention to the 
atmospheric factors affecting dry weather (at least for heat waves), e.g., the cyclones and anticyclones 
that control cloud cover over most of the northern extratropics. 

b) Persistent anticyclones 
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The key role of anticyclones in generating Eurasian drought and heat waves has important 
implications for their spatial structures and time scales. In Eurasia, severe drought conditions in one 
region (say, European Russia) are at times accompanied by wet and cool conditions to the west over 
Europe and/or to the east over parts of Siberia. This is suggestive of an east- west wave structure 
underlying the surface temperature and precipitation anomalies and thus of strong atmospheric 
controls. Such structures are evident in the analysis of Eurasian heat waves provided by Gershunov 
and Douville (2008). They note that “In both model and observations, there is a strong interannual 
propensity for Far Eastern Europe to be cold during heat wave summers in West-Central Europe. Both 
recent extreme European heat wave summers of 1994 and 2003 were cold in far-eastem Europe and 
warm over north-central Siberia, thus exhibiting Eurasian summer temperature wave train conditions 
typical of large European heat waves.” Stankunavicius et al. (2012) carried out an empirical 
orthogonal function (EOF) analysis of surface air temperature (SAT) and sea level pressure (SLP) for 
every (2 month) season over Eurasia based on NCEP/NCAR reanalyses for the second half of the 20 th 
century. They found clear evidence of wave structures in the leading modes of SAT variability during 
early and late summer. Sato and Takahashi (2006) identified a southern Eurasian wave train 
extending far enough eastward to affect Japan. 

Using NCEP/NCAR reanalysis data, Stefanon et al. (2012), with a clustering approach, identified six 
major types of European heat waves for the period 1950-2009. The types found include a Russian 
cluster, a Scandinavian cluster, a western European cluster, and an eastern European cluster, with the 
temperature anomalies in phase with the anticyclonic (positive 500mb geopotential height) anomalies. 
They found that drought appears to be a pre-requisite to heat wave occurrence in western and eastern 
European heat waves (rainfall deficits in southern Europe), but not for the more northerly Russian or 
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Scandinavian heat waves (see Section 3 below). 


c. The leading modes of surface temperature and precipitation covariability 
The above studies on drought and heat wave characteristics suggest that we can effectively quantify 
hydro-dynamical variability over Eurasia in terms of the combined monthly temperature and 
precipitation variability. In order to do this efficiently, we employ a rotated 8 empirical orthogonal 
functions (REOFs) analysis. The basic quantities used in the calculation of the REOFs (the 
normalized monthly temperature and precipitation fields) are the same as those used in the calculation 
of the Ped (1975) drought index (Section 2a). We focus here on monthly rather than seasonal means 
to better capture the variability associated with persistent large-scale atmospheric waves. 

The first REOF (top panels of Fig 3) shows a clear wave structure in both the temperature and 
precipitation, with anomalies of alternating sign spanning Eurasia (both across the north over Siberia 
and to the south into China). The greatest temperature loading (shown here as positive) is centered on 
European Russia (the European Plain) west of the Ural Mountains. The associated precipitation 
loadings have negative values on the southeastern quadrant of the main wann anomaly (just north of 
the Caspian Sea), suggestive of a dynamical link between the temperature and precipitation anomalies. 
Positive values for precipitation occur over central Europe, Scandinavia, northern Siberia, and 
mountains of Central Asia. The corresponding time series of the leading REOF (referred to as the 
rotated principal component, or RPC) in Figure 4 show that this pattern is associated with a trend 
towards more positive values over the last thirty years; it thus appears to have played an important 
role in the Russian heat wave of 2010 (relatively large positive values for June, July and August). 

8 Rotation (Richman 1986) acts to spatially localize anomalies, and has been found by the authors to produce more physically 
realistic patterns of variability compared with unrotated EOFs. We note that the REOF methodology has no inherent tendency to 
produce wave structures, in fact, the localization would tend to de-emphasize connections at large distances. 
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The second REOF again shows a wave structure, with the largest loading in the temperature field just 
east of the Urals, indicating an east- west phase shift with respect to the first REOF. Positive values 
also occur over southern Europe, while negative values occur over Scandinavia, northeastern Europe, 
and much of eastern Asia. The associated precipitation loadings show negative precipitation 
anomalies just to the east of the Ural Mountains (again on the southeastern quadrant of the main 
positive temperature anomaly) and over much of Europe, while positive anomalies stretch from 
Scandinavia southeastward across Kazakhstan, Mongolia and China. The associated RPCs show that 
this pattern was very pronounced in August of 2003, during the height of the 2003 European heat 
wave. We note that the first REOF was also very pronounced (negative) in 2003 though this occurred 
in June of that year at the start of the heat wave. The associated RPC shows no clear trend in the last 
three decades (Fig. 4). 

REOFS 3-5 are also characterized by wave structures, with the maximum positive temperature 
anomalies centered just east of the Caspian Sea, just north of Mongolia, and over northern Europe, 
respectively. REOF 3 differs somewhat from the others in that it is indicative of a more southerly 
wave path. In all three cases, the main negative precipitation anomalies are either in phase or slightly 
to the east/southeast of the main positive temperature anomalies. The associated RPCs (particularly 
RPC3 and RPC4) indicate a change towards more positive values after about 1995. We will come 
back to the trends in Section 4. In the following section, we focus on the mechanisms responsible for 
such wave structures. 

3. Physical Mechanisms 
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a. A Review 


As mentioned in the previous section, it has long been recognized that persistent anticyclones play a 
fundamental role in the generation of drought and heat waves over northern Eurasia. Buchinsky 
(1976) summarizes some of the key aspects of droughts in the central part of European Russia and the 
Volga region, noting that over 70% of them are associated with persistent anticyclones that act to 
disrupt the predominantly zonal flow and eastward progression of weather systems. He notes that 
these are primarily Arctic anticyclones that advance from the Barents or even the Kara Sea and 
become stationary over the plains. Similarly, the work of Selianinov (1928) and others, as 
summarized in Kleshenko et al. (2005), showed that drought in the arid regions of Russia and other 
Commonwealth of Independent States (CIS) regions results from the penetration of anticyclonic air 
masses from the Artie. They note that these anticyclones can act in concert with anticyclones at the 
southern (40-50°N) and high (approx. 75°N) latitudes. The former become more important further to 
the west where, for example, Ukraine is impacted by the Azores high. They note that most often the 
Arctic and Azores intrusions are combined in the Lower Volga and the Southern Yuzhny Ural, leading 
to pronounced drought conditions. 

Similarly, eastern Europe (e.g., Poland) periodically experiences drought related to a persistent 
stationary anticyclone (an east European high) that joins with the Azores anticyclone via central 
Europe (Farat et al. 1998). As noted by Golubev and Dronin (2004) “An especially strong drought 
takes place when an anticyclone is fed by an air mass from an Azores anticyclone moving in from the 
West. Moving across Europe, the air mass loses its humidity and reaches European Russia completely 
dry (Protserov, 1950). The droughts resulting from these large scale atmospheric processes usually 
cover vast territories of Russia, including the Northern Caucasus, the Middle and Low Volga basin, 
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the Urals, and periodically spread over the central chernozem region and even the northern regions of 
European Russia. For example, the drought of 1946 covered 50 percent of total agricultural land of the 
USSR. As a result, the scale and consequences of droughts can be catastrophic for the country.” 

The physical mechanisms that determine the persistence and scale of the northern Eurasian 
anticyclones are still not well understood, though atmospheric blocking has long been considered 
important. Studies of blocking that focus on the Atlantic and impacts on Eurasia go back to Obukhov 
et al. (1984) and a number of earlier studies reviewed therein. That study in particular reviewed 
various potential mechanisms of blocking, including those linked to orography and the instability of 
the Polar Jet, and it emphasized atmospheric blocking as a precondition for drought in summer, with 
both the downward movement of air within the associated anticyclone (acting to heat and dry the air) 
and the blocking of the westerlies (inhibiting the inflow of moisture from the west) contributing to the 
drought conditions. More recently, Nakamura et al. (1997) contrasted Pacific and Atlantic blocking 
events and found that incoming wave activity associated with a quasi- stationary Rossby wave train is 
of primary importance in the development of blocking over Europe, while the forcing from synoptic 
scale transients is key for the development over the North Pacific. 

In addition to blocking, a number of other large scale modes of variability can affect northern Eurasia 
on weekly to monthly time scales. The important role of the Northern Annular Mode (NAM) for 
Eurasian climate has been documented in numerous studies (e.g., Thompson and Wallace 2001). 

While many studies have focused on the winter season, others have documented the impact of the 
NAM on variations in land surface phenology such as the start of the growing season and the timing 
of the peak NDVI over northeastern Russia (e.g., De Beurs and Henebry 2008). Rocheva (2012) 
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linked the persistent 500mb height anomalies over European Russia and western Siberia to the East 
Atlantic/Western Russia (EA/WR) and Scandinavian (SCA) patterns of variability, respectively, from 
May through July - patterns that are still not well understood. 

Bothe et al. (2003) li nk ed drought over Tibet to a Eurasian wave train that spans Eurasia from 
Scandinavia to the South China Sea. They associated the development of the wave to strong 
anticyclonic activity over northern Europe/Scandinavia, which in turn is supported by anomalous 
transient eddy activity associated with the North Atlantic stonn track. Ding and Wang (2005) 
identified a wave number 5 summertime circumglobal teleconnection pattern confined to the summer 
jet wave guide with significant impacts on interannual (and intraseasonal - Ding and Wang 2007) 
temperature and precipitation variations over much of Eurasia and North America, apparently 
maintained by heat sources associated with the Indian monsoon. Schubert et al. (2011) examined the 
role of stationary Rossby waves on intraseasonal summertime variability in the Northern Hemisphere 
extratropics and showed that many of the extreme events, including the 2003 European and 2010 
Russian heat waves, are associated with a particular recurring Eurasian stationary wave pattern that 
affects much of the northern Eurasian continent. This, along with other summertime wave structures, 
were found to be primarily forced by sub-monthly vorticity transients, though it was also found that 
the waves do at times contribute substantially to the seasonal mean anomalies, suggesting some 
impact from other longer-term (e.g. SST) forcing. 

Uncertainties about the causes of persistent northern Eurasian anticyclones result from limitations in our 
understanding of the basic dynamical mechanisms involved and from uncertainties about the impact of 
global warming, especially in regard to the occurrence of some of the most extreme events (e.g., Dole et al. 
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2011; Schneidereit et al. 2011; Lau et al. 2012; Galarneau et al. 2012; Trenberth and Fasullo 2012; Lupo et 
al. 2012). On the one hand, for example, Dole et al. (2011) emphasized the important role of internal middle 
latitude atmospheric dynamics in producing an intense and long-lived blocking event and associated 
anticyclone (producing the wannest July since at least 1880 in western Russia), and they concluded that 
neither human influences nor slowly varying ocean boundary conditions contributed substantially to the 
magnitude of the event. They also stated that “severe drought occurred with the Russian heat wave, making 
it likely that land surface feedbacks amplified this heat wave’s intensity.” Trenberth and Fasullo (2012), in 
contrast, linked the unusual anticyclone to the development of a large-scale Rossby wave train - suggesting 
that the wave train was forced by anomalous convection in the tropical Atlantic and northern Indian oceans. 
They also argue that the heat wave intensified through the cumulative impact of local land feedbacks, linked 
to increased greenhouse gases. Lau and Kim (2012) highlighted the role of this wave train in linking the 
Russian heat wave to the Pakistani floods, with land feedbacks acting to amplify the Russian heat wave, and 
moisture transport from the Bay of Bengal (associated with the northward propagation of the monsoonal 
intra-seasonal oscillation) helping to sustain and amplify the Pakistani rains. They argue that the western 
Russian blocking event was itself instrumental in forcing the Rossby wave. Galereau et al. (2012) 
highlighted the importance of circulation around the blocking ridge accompanied by enhanced subsidence in 
the intensification of the heat wave. They also found that downstream energy dispersion from source regions 
over the North Atlantic modulated the structure and intensity of the blocking anticyclone over western 
Russia. 

Schneidereit et al. (2012) argue that a number of factors at several different time scales were at work 
during the 2010 heat wave. They show that the shift to La Nina conditions modulated the stationary 
wave pattern, supporting the blocking high over Eastern Europe. Also, they found that a polar Arctic 
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dipole mode projected on the mean blocking high, and that transients acted to maintain it. At 10-60- 
day time scales they identified three different paths of wave action that also contributed to the 
persistent blocking conditions. 

While numerous studies have addressed the important role of soil moisture feedbacks in European 
droughts (e.g., Ferranti and Viterbo 2006; Seneviratne et al. 2006; Fischer etal. 2007; Vautard et al. 
2007; Zampieri et al. 2009; Stefanon et al. 2012), far fewer studies have focused on soil moisture 
impacts in the rest of Eurasia. Cherenkova (2012) examined the precursors to summer drought in the 
European Territory of Russia, finding that of the five most extensive hazardous droughts that occurred 
between 1936-2010 (1936, 1938, 1972, 1981 and 2010), three (1936,1938 and 1972) were preceded 
by dry winters and springs, which created conditions for further drought development in the summer. 
The 1981 drought was not preceded by a dry winter and spring, and they suggest that this explains the 
smaller area covered by that drought. The 2010 drought was preceded by a cold winter without 
precipitation deficits, but they suggested that the cold temperatures did impact the snowmelt and 
spring soil moisture deficits. As already mentioned, Lau and Kim (2011) found that the 2010 Russian 
heat wave was amplified by the underlying extensive region of dry soil conditions. 

Hirschi et al. (2011) examined the relationship between soil moisture, drought, and summer heat for 
central and southeastern Europe, based on observational indices for 275 station observations. They 
found that dry soil conditions intensified hot extremes in the southeastern (Romania and Bulgaria) 
area, especially for the high end of the distribution of temperature extremes, whereas this was not the 
case for central Europe (Austria and Czech Republic); they further noted while that the fonner area is 
characterized by soil moisture-limited evaporation, the latter is characterized by energy-limited 
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evaporation. Mueller and Seneviratne (2012) show that the dryness-temperature relationship is 
important in many areas of the world including much of Eastern Europe (extending east into European 
Russia to about 50°E), where the probability of occurrence of an above-average number of hot days 
with preceding precipitation deficits is over 60%. Volodin (2011) analyzed the causes of “super- 
extreme” anomalies of summer surface air temperature in a suite of GCM experiments and reanalyses, 
focusing on the summer 2010 hot spell over ER as well as similar hot spells in Western Europe (2003) 
and the contiguous United States (in 1980 and 2007). He showed that, in addition to the atmospheric 
factors acting during the peak month of drought (in the case of ER in July 2010 this was a prolonged 
atmospheric blocking event), preceding monthly anomalies of soil moisture located windward of the 
drought significantly enhanced the temperature anomaly. This behaviour repeated itself over ER in 
the summer of 2012 where the drought (and the soil moisture anomaly) began initially over 
Kazakhstan and the southernmost areas of ER and gradually expanded northward. Lorenz et al. 

(2010) analyzed regional climate model simulations to show that soil moisture memory also acts to 
increase the persistence (in addition to the intensity) of heat wave events. 

Koster et al. (2006b, see their Figure 11) showed that the observed spatial pattern of interannual JJA 
temperature variance over North America can be reproduced by an AGCM only when soil moisture 
feedback processes are allowed to operate in the model, a strong indication that soil moisture 
variability contributes significantly to temperature variability. An analogous figure for Eurasia is 
shown here in Figure 5. Figure 5a, from observations (MERRA), shows monthly temperature 
variance for the JJA period over 1980-2012. To first order, the free-running GEOS-5 AGCM (Figure 
5b) reproduces this structure, with high variability in the most northern parts of Eurasia and another 
band of high variability centered at about 50°N. (The variances produced by the free-running model 
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and MERRA differ mostly in their amplitudes, with weaker values seen in the fonner.) Figure 5c 
shows the temperature variances generated by the GEOS-5 AGCM when soil moisture feedback 
processes are artificially disabled, a condition achieved here by continually resetting the land model’s 
soil moisture prognostic variables to seasonally-varying climatological values (see Appendix A). 
Disabling soil moisture feedback significantly reduces the variances along a swath through the center 
of the continent, extending from southern Europe eastward across the Caucasus to Kazakhstan, 
Mongolia, and northern China (Figure 5d). This swath of reduction is indeed where we expect it to 
be, located at the transition between the wet climate to the north and the dry climate to the south; 
evaporation variance associated with soil moisture variations tends to be maximized in such a 
transition regime (see Koster et al. 2006b). In the AGCM, soil moisture feedback is unequivocally 
responsible for enhanced temperature variance along this swath, and we can speculate that the same is 
true in nature. We will come back to the role of the land later in our discussion of long-term trends 
and predictability. 

The role of SST in seasonal to decadal climate variability over Eurasia is also still not well 
understood. Again, much of the analysis of the role of SST has focused on impacts in Europe, though 
a number of these studies have implications for regions to the east. Ionita et al. (2012) analyzed the 
self-calibrating Palmer drought index (van der Schrier et al. 2006) for the period 1901 - 2002 and 
found considerable interannual and multi-decadal variability in summer moisture over Europe that 
was tied to SST variability. In addition to a drying trend over Europe associated with warming SST 
over all oceans, they found a link between previous winter La Nina and negative PDO events and 
summer dry conditions over southern Europe extending into western Russia, with wet conditions over 
the Scandinavian Peninsula, with the atmospheric anomalies resembling aspects of the PNA and (the 
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positive phase of the) NAO. They also found a li nk with the cold/negative phase of the AMO that 
leads to summer drought conditions extending across southern Scandinavia, southeastern Europe and 
into northwestern Russia. They cite the extremely hot and dry summers of 1921 and 1972 over the 
central and northern regions of Russia (Buchinsky 1976) as examples of events that coincide with an 
AMO in its negative phase. 

Sedlacek et al. (2011) hypothesized that the SST anomalies in the Barents and the Arabian Seas 
combined to produce warming over Eurasia during 2010, thus contributing to the heat wave; they 
suggest that such a dynamic response to SST (in particular to the expected warming and reduction in 
sea ice over the Barents Sea) will contribute to more frequent heatwaves over Eurasia in the future. 
Wu et al. (2012) examined the impact of the North Atlantic Oscillation on the relationship between 
the East Asian summer monsoon and ENSO and found, among other things, that an anomalous spring 
NAO induces a tripole SST anomaly in the North Atlantic which persists into summer and excites 
downstream development of a Rossby wave train that modulates the blocking highs over the Ural 
Mountains and O kh otsk Sea. While the main impact of Arctic Sea ice reduction occurs during winter 
(Deser et al. 2007), a recent observational study by Francis and Vavrus (2012) suggests that the 
reduction in Arctic Sea ice slows the progression of Rossby waves by weakening the zonal winds and 
increasing wave amplitude. They argue that while these impacts are strongest during winter and 
autumn, they are also apparent in summer (possibly due to earlier snow melt on high-latitude land) 
and therefore contribute to more extreme summer weather events including Eurasian heat waves. 

b. The role of Stationary Rossby waves 
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A recurring theme in the above discussion is the role of Rossby waves. Ambrizzi et al (1995) provide 
one of the first studies to isolate the teleconnectivity associated with the boreal summer wave-guides 
and preferred wave propagation patterns towards and away from the wave-guides. Again, Schubert et 
al. (2011) identified a particular recurring Rossby wave (forced by submonthly vorticity transients) 
that extends across northern Eurasia and that contributes significantly to monthly surface temperature 
and precipitation variability, playing an important role in the generation of the 2003 European and 
2010 Russian heat waves. 

To the extent that Rossby waves are an important component of summer Eurasian temperature and 
precipitation variability, we would expect that the leading surface temperature and precipitation 
REOFs shown in Figure 3 would be tied to such atmospheric waves. The correlations between the 
leading RPCs and the monthly 250mb v-wind (Figure 6) suggest that this is indeed the case. The 
correlations with the first two RPCs show two clear wave structures that are approximately in 
quadrature extending across northern Eurasia. In fact these closely resemble the Schubert et al. (2011) 
basic wave structure of the leading REOF of the monthly 250mb v-wind mentioned above. The 
correlation pattern associated with the first RPC differs somewhat from that of the second in that the 
anomalies seem to extend around the globe, and there is a clear signature of a split in the wave over 
Europe with the northern component extending across Eurasia to the north of the mean jet, and a 
southern component that appears to use the mean jet as a wave guide (this correlation pattern is very 
similar to the actual 250mb v-wind anomalies during July 2010 - cf. Fig. 7). The correlations 
associated with the second RPC suggest a wave development that is more confined to the northern 
part of Eurasia (north of the mean jet), and resembles the June 2003 v-wind anomalies (Fig. 7). 
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The correlations with the third RPC (Figure 6) show a wave structure that is more confined to the 
mean jet throughout the Northern Hemisphere; over Eurasia it appears to affect primarily southern 
Europe and the regions east of the Caspian Sea. This pattern dominates the v-wind anomalies, for 
example, in August 1992 (Fig. 7). The fourth and fifth RPCs are associated with wave structures 
similar to those of the other leading modes, but with that of the fourth having its largest amplitude 
over the eastern half of Eurasia (e.g., August 2001 in Fig. 7), and that of the fifth having its largest 
amplitude over northeast Atlantic and northern Europe (e.g., July 1994 in Fig. 7). 

The potential role of SST anomalies in forcing the leading REOFs was examined by computing the 
simultaneous and time-lagged correlations with the global monthly SST anomalies (Figures not 
shown). The results indicate that the correlations are generally weak (absolute values less than 0.3). 
An exception to that is RPC1, which has somewhat larger negative correlations (between -0.3 and - 
0.4) in the tropical eastern Pacific at both 0 and -1 month lags, suggesting a weak link to ENSO. Also, 
RPC 3 has positive correlations with SST (between 0.3 and 0.4 at lag 0) over the North Atlantic, with 
a similar pattern of correlations (but weaker) occurring at -1 lag. The largest correlations with SST 
occur for RPC5 (values greater than 0.5 at lag 0) over the far eastern North Atlantic and 
Mediterranean Sea in the immediate vicinity of Europe: these likely reflect the response of the SST to 
the changes in atmospheric forcing associated with the wave itself. The above results indicate that 
SSTs have only a weak (if any) impact on the development of these waves on monthly time scales, 
with perhaps ENSO and the North Atlantic SST having some influence on RPCs 1 and 3, respectively. 
An important caveat here is that the above correlations reflect primarily interannual linkages in the 
monthly statistics, rather than subseasonal linkages. In fact, if we remove the interannual component 
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of the variability, the correlations with SST are even weaker for all RPCs except for the simultaneous 
correlations associated with RPC 5 in the vicinity of Europe. 

The forcing of such waves by submonthly vorticity transients is illustrated in Fig. 8 in the context of a 
stationary wave model (SWM- Ting and Yu 1998) forced by an idealized localized vorticity source in 
the North Atlantic jet exit region (see Schubert et al. 201 1 for details). The results show that an 
atmospheric wave structure very similar to that associated with REOF 1 develops in the SWM within 
about three weeks. While such Rossby waves (driven by internal atmospheric dynamics) appear to be 
a ubiquitous component of summertime weekly to monthly atmospheric variability over Eurasia, the 
mechanism that leads to their occasional persistence and amplification are as yet unclear. An 
assessment of the potential role of soil moisture and a further assessment of SST forcing will be made 
in the following two sections, where we examine longer-term (seasonal to decadal) variations and the 
predictability of such extreme events. 

4. Past Long term Behavior and Trends 

a) A Review 

The 2003 European and 2010 Russian heat waves, in addition to prompting numerous papers on 
causes and impacts, highlighted the on-going debate about whether such events are early 
manifestations of global warming. For example, Rahmstorf and Coumou (2011), employing a 
stochastic model to examine the effect of warming trends on heat records, concluded that, with a 
probability of 80%, “the 2010 July heat record would not have occurred” without the large-scale 
climate warming seen since 1980, most of which has been attributed to the anthropogenic increase in 
greenhouse gas concentrations. In contrast, as already mentioned, Dole et al. (2011) conclude from 
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their analysis of dynamical mechanisms that neither human influences nor slowly varying ocean 
boundary conditions contributed substantially to the magnitude of the 2010 event. Otto et al (2012) 
examined the results from a large ensemble of atmospheric general circulation model simulations and 
concluded that “ there is no substantive contradiction between these two papers, in that the same event 
can be both mostly internally-generated in terms of magnitude and mostly externally-driven in terms 
of occurrence-probability. The difference in conclusion between these two papers illustrates the 
importance of specifying precisely what question is being asked in addressing the issue of attribution 
of individual weather events to external drivers of climate.” 

In addition to the current debate on whether the nature of extreme events is changing, there is also 
ongoing debate about basic trends in both the mean precipitation and surface temperature. 
Meshcherskaya and Blazhevich (1997) used station data to study changes in drought over the 
European and Asian parts of the former Soviet Union (FSU) for the period 1891 - 1995. They found 
that trends in their drought and excessive moisture index (DM, for May- July) are statistically 
significant only in the Asian part of the FSU and that the increased dryness is largely the result of 
temperature increases, with a small but statistically significant contribution coming from a decrease in 
precipitation. These results have recently been updated and expanded by Groisman et al. (2013), who 
showed that while heavy rainfall frequencies have increased in the past two decades, mean 
precipitation has grown more slowly or has even decreased, with an accompanying increase in the 
frequency of no-rain periods over most of Northern Eurasia south of 60°N. 

Alexander et al (2006), using updated station data from the more recent record (1951-2003), examined 
a number of climate indices (see also Frich et al. 2002) and found significant changes in extremes 
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associated with warming. In particular, they found that much of Eurasia is characterized by a 
significant decrease in the annual number of cold nights and an increase in the number of warm 
nights. These results hold for all seasons, with the largest changes occurring during MAM and the 
smallest during SON. Corresponding behavior is also seen in a subset of stations with records going 
back to 1901. Precipitation indices show a tendency toward wetter conditions throughout the 20 th 
century. 

Frey and Smith (2003) examined precipitation and temperature trends in station observations from 
western Siberia, a region with a large percentage of the world’s peatlands, and one that contributes 
substantially to the terrestrial freshwater flux into the Arctic Sea. They found robust patterns of 
springtime wanning and wintertime precipitation increases, with the AO playing an important role in 
non-summer wanning trends. As noted by the IPCC (2001), the AO (and NAO) had been in phase 
since the 1970s, producing enhanced westerlies and extratropical cold season wanning across much of 
Eurasia. 

Batima et al. (2005) examined data from 60 meteo-stations spanning Mongolia for the period 1940- 
2001 and found that the mean annual surface temperature has risen by 1.66°C over the 62-year period, 
warming faster in winter than summer. The warming is more pronounced in mountainous areas and 
their valleys and is less pronounced in the Gobi desert. They also find a statistically insignificant 
decrease in annual mean precipitation, with winter and spring showing a decrease but summer and fall 
showing no change. Even without clear evidence for an increase in summer temperatures, summer 
heat wave duration has increased by 8-18 days, depending on location. 1998 was the warmest year of 
the last century, and Mongolia experienced drought for the next four years (1999-2002). Batima et al. 
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(2005) further note that the intense drought spells in recent years are most likely the result of both 
increased temperature and decreased precipitation. They emphasize that the environment and climate 
plays a key role in the sustainability of Mongolia - animal husbandry employs 47.9% of the total 
population, producing 34.6 percent of the agricultural gross production, and accounting for 30% of the 
country’s exports. Nandintsetseg et al. (2007) found an almost 2°C increase in temperature in 
northern Mongolia between 1963-2002, along with a significant increase in wann extremes and a 
decrease in cold extremes. On average, they found neither a significant decrease in the maximum 
number of consecutive dry days nor an increase in the number of wet days 

Robock et al. 2005 examined 45 years (1958-2002) of soil moisture observations over Ukraine and 
found an increase in soil moisture over those years, despite a slight warming and a decrease in 
precipitation. They suggested that this is the result of increased aerosols in the troposphere leading to 
decreased solar insolation, which acts to reduce evaporation; the reduced evaporation in turn leads to 
increased surface temperature and soil moisture. 

IPCC (2007b, Table 10.2) summarize some of the key trends in northern Eurasia, with Russia 
experiencing a 2 to 3°C rise in the past 90 years that is most pronounced in spring and winter. 

Changes in precipitation in Russia are highly variable with a decrease during 1951 to 1995, and an 
increase in the last decade. Central Asia experienced a 1-2 °C rise in temperature per century, with no 
clear trend in precipitation between 1900 and 1996. Mongolia has seen a 1.8°C increase in the last 60 
years that is most pronounced in winter; Mongolian precipitation has decreased by 7.5% in summer 
and has increased by 9% in winter. 
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Analyses covering longer time periods are also available. Briffa et al. (1995) report on a 1000-year 
tree-ring reconstruction of summer temperatures over the northern Urals; they show that the mean 
temperature of the 20 th century is higher than that of any other century since AD 914. Demezhko and 
Golovanova (2007) reconstructed ground surface temperatures from AD 800 onward based on 
borehole temperature logs and 170 years of meteorological data over the southern and eastern Urals. 

They conclude that the mean temperature during the Medieval maximum (1 100-1200AD) was 0.4°K 
higher than that for the period 1900-1960. They also conclude that cooling during the “Little Ice Age” 
culminated in about AD 1720 with a mean surface temperature 1.6°K below the 1900-1960 mean, and 
they note that the contemporary warming began about a century prior to the first instrumental records 
in the Urals, with the mean rate of warming increasing in the final decades of the 20 th century. 

The recent Special Report of IPCC on extremes (IPCC, 2012) provides an updated summary of the current 
confidence placed in recent trends of heat waves and droughts. The report notes that in Asia there is “overall 
low confidence in trends in dryness both at the continental and regional scale, mostly due to spatially varying 
trends, except in East Asia where a range of studies, based on different indices, show increasing dryness in 
the second half of the 20th century, leading to medium confidence.” They also note that since 1950, there is 
medium confidence in a wanning trend in daily temperature extremes over much of Asia. 

b) A model-based analysis of recent trends (1979-2012) 

In this section, we utilize numerical simulations to provide further insight into the nature of recent 
variability and trends over Eurasia. These simulations take the form of full global reanalyses 
(MERRA and ERA Interim), AMIP-style simulations using the GMAO GEOS-5 system, and 
simulations with more idealized SST forcing. 
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One of the intriguing aspects of RPC 1 in Fig. 4 is the apparent trend or shift in the time series from 
being predominantly negative prior to about 1995 to predominantly positive thereafter. There is an 
indication of a similar shift in RPCs 3 and 4. The first three panels on the left of Fig. 9 (derived from 
the two reanalyses and from station observations) indicate that these shifts appear to be part of a 
hemispheric-wide pattern of warming over the last three decades, with the maxima over Eurasia 
centered over European Russia and Mongolia/eastern Siberia. While the maps from the reanalyses 
differ somewhat from that constructed with the station observations, especially regarding the 
amplitude of the changes, overall they agree on the main regions of warming. As for precipitation 
(the first three panels on the right side of Fig. 9), the patterns of change are more complex, with 
decreases covering parts of northeastern Europe, European Russia, Kazakhstan, southeastern Siberia, 
Mongolia, and northern China, and with increases found across Siberia north of about 60°N. 

The extent to which the above trends are a reflection of global warming and/or the result of other 
long-tenn (decadal) variability is still an open question. Some insight into this issue can be gained 
from the analysis of free-running climate model simulations. We examine now an ensemble of 12 
GEOS-5 AMIP simulations driven by observed SST and GHG forcing over the period 1871 -present. 
The 1996-2011 minus 1980-1995 differences for the ensemble mean are shown in the bottom panels 
of Figure 9. Overall the model results are consistent with the reanalyses and observations, showing 
warming over basically the same regions across Eurasia (southern Europe and European Russia, 
Kazakhstan/southem Siberia, Mongolia, and northern China) though with weaker amplitude. We 
note, however, that individual ensemble members (not shown) exhibit changes as large as the 
observed and that there is substantial intra-ensemble variability in the detailed spatial patterns of the 
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differences, with some showing the same two-lobed structure found in the observations. The model 
ensemble mean also reproduces to some degree the overall pattern of precipitation changes (though 
again with weaker amplitude than observed), including the tendency for precipitation deficits over 
European Russia and over Siberia south of about 60°N, and for precipitation increases to the north. 

The AMIP results in Figure 9 suggest that SST variations and perhaps the direct GHG forcing are 
contributing significantly to the observed JJA trends in Eurasian surface temperature and precipitation 
seen over the last three decades. Figure 10 (left panel) shows the linear trend in observed SST during 
that period (1980-2011). The SST trend pattern shows aspects of overall warming combined with a 
La Nina-like pattern in the Pacific and a positive Atlantic Multi-decadal Oscillation (AMO) pattern in 
the Atlantic. This is compared (right panel of Fig. 10) with one of the idealized SST forcing patterns 
used recently by the USCLIVAR Drought Working group to force several different climate models 
(Schubert et al. 2009). This latter pattern is the sum of the three leading REOFs of annual mean SST, 
consisting of a PDO/La Nina -like pattern, an AMO-like pattern and the warming trend pattern. With 
the exception of the Indian Ocean, the similarity of this idealized pattern to the recent (three decade 
long) trend pattern is striking, suggesting that the recent trends are a mixture of both decadal 
variability and long-term trends. 

Figure 1 1 shows, for the average of three of the models that participated in the USCLIVAR Drought 
Working group project (Schubert et al., 2009) and GEOS-5, the JJA surface temperature response to the 
idealized SST pattern shown in the right panel of Fig. 10. The results are based on 50-year long simulations 
for all the models except CFS which was integrated for 35-years. The models produce warming (top left 
panel) over most of Eurasia between 30-60°N. The precipitation anomalies (top right) consist of deficits 
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over central Eurasia (centered on about 50°N, 95°E) and parts of Europe. Positive precipitation anomalies 
occur over much of the northern regions of Russia, especially east of about 70°E, extending into northeastern 
China. Additional runs with these models (not shown) indicate that the Pacific and Atlantic SST patterns act 
to focus the warming and precipitation deficits in the middle latitude band between 30-60°N, as well as to 
produce some regional (east/west) variations that differ from model to model. The SST trend pattern acts to 
expand and enhance the regions of warming, with an overall tendency to warm the continents everywhere. 9 
These results are generally consistent with those shown in Fig. 9, supporting the idea that the main features 
of the northern Eurasian precipitation and temperature trends of the last three decades are largely forced by 
the leading patterns of SST variability (the global trend and the two dominant patterns of SST variability in 
the Pacific and Atlantic Oceans). The bottom panel of Figure 1 1 suggests that the surface warming and 
precipitation changes are linked to a tendency for all the models to produce a band of positive upper 
tropospheric height anomalies throughout the middle latitudes of both hemispheres in response to the 
imposed SST patterns, and that (as shown by additional runs isolating the SST trend impacts - not shown) 
these positive height anomalies, while basically forced by the Pacific and Atlantic SST patterns, are 
amplified with the additional forcing of the SST trend pattern. 

c) An Analysis of Long-Term Variability (1871-2010) 

To put the trends of the last three decades in perspective, we now turn to temperature records going 
back to the late 19 th century. The left panels in Figure 12 show the time series of JJA mean 
temperature for the period 1871 to 2010, based on the CRUTEM4 data, for four different regions in 
northern Eurasia: Europe, European Russia, south-central Siberia, and a cold desert region just east of 
the Caspian Sea, centered on the Aral Sea. (The definition of the regions was guided by the regions of 

9 There are some regions of cooling over land in response to the SST trend pattern (which itself has some spatial variability) 
though these tend to be relatively small in area and have small amplitude of cooling. 
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maximum T2m loadings of the leading REOFs shown in Fig. 3.) All four regions show 
predominantly positive anomalies beginning shortly after 1990, though this is most pronounced for the 
European region. The 2010 heat wave stands out in the European Russia time series, though there are 
some other years with large anomalies, including 1972 (during the “100-year” drought; see also 
Appendix B). While there is substantial interannual and decadal scale variability in all of the time 
series, there is also evidence of a long-tenn positive trend, though the trend values appear to depend 
somewhat on the observations during the late 19 th century, which are likely not very reliable in the 
CRUTEM4 data set. This is illustrated through comparison with another data set (NOAA MLOST - 
right panels of Fig. 12), indicating some differences during the early years, especially for the more 
eastern region. The latter dataset shows clearer long-term trends in part because it does not include 
the 1870s, which in the CRUTEM4 data is a period of positive anomalies. Note that extensive 
standard surface air temperature observations over the Russian Empire territory began in 1881 
(Vannari 1911). 

Figure 13 is the same as Figure 12, but constructed from the output of two representative members of 
the aforementioned 12-member ensemble of 130-yr AMIP simulations with GEOS-5 (see Appendix 
A). Each time series shows a basic character that is remarkably similar to that of the observations, 
with a shift towards positive anomalies starting in the 1990s. The long-tenn trend in all four regions 
is, in fact, even more pronounced in the model simulations. The simulations also show a few very 
extreme anomalies. In particular, we point out the unusually large positive (+3°C) anomaly simulated 
in 2001 in European Russia in one of the ensemble members (second row of Fig. 13, on left). This 
event has a temperature signal comparable in magnitude to that observed during the 2010 Russian heat 
wave (Fig. 12) and will be examined in more detail in the next section. We note that similar events 
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(temperature anomalies near +3°C) occur in some of the other ensemble members, though in general 
they are quite rare and are limited to the recent decades (e.g., there were only three such events in 
Europe and only two in European Russia in the entire ensemble of 12 runs, i.e., over a span of 1680 
simulation years). 

We next examine whether the interannual JJA mean surface temperature (T s ) variations in these 
regions are li nk ed to SST variability. Figure 14 shows the temporal correlations of the regional mean 
T s for the European Russia region with T s values everywhere across the globe (using SST over the 
oceans). The calculation is limited to the years 1901-1980 in order to avoid the earliest years with 
little observational data and to avoid the three most recent decades, which show the shift towards 
positive values. In addition, a linear trend was removed from all time series prior to computing the 
correlations. The results for the observations (top panel) 10 show a wave structure over northern 
Eurasia that is very similar to the leading REOF of the monthly data for the recent three decades (Fig. 

3 top left panel), suggesting that the seasonal data also project onto a Rossby wave structure, 
consistent with the findings of Schubert et al. (2011). There are also positive correlations over the 
Atlantic and the eastern Pacific, suggesting some link to the SST. The results for the model are shown 
in the middle and bottom panels of Figure 14. The middle panel shows the correlations computed 
separately for each ensemble member and then averaged over the 12 ensemble members - a result that 
is more comparable to the correlations based on the observations. The model results clearly show the 
same basic wave structure of the correlations over northern Eurasia. By this measure, the link to the 
oceans is weak, with only small positive correlations (0. 1-0.2) that are mainly confined to the Atlantic, 
however, there are individual ensemble members (not shown) that have correlations resembling those 
based on the observations (highlighting a considerable unforced component to the observational 
10 Absolute values greater than 0.22 are significant at the 0.05% level ( www.mtsu.edu/~dwalsh/436/CORRSIG.pdf) 
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results). The bottom panel in Figure 14 shows the correlations with the ensemble mean T s . This 
calculation isolates the impact of the forcing common to all the ensemble members (SST and GHGs), 
showing, for example, correlations over the North Atlantic that exceed 0.4. Also, it is noteworthy 
that the pattern of correlations over the eastern Pacific is very similar to that based on the 
observations. Other regions with relatively large correlations include Northern Africa, southern 
Eurasia, Canada and the western US. Further analysis (not shown) indicates that the patterns of the 
correlations are quite sensitive to the location of the target area. Analogous global maps of T s 
correlations with T s values in southern Europe, for example, show substantial negative correlations in 
the tropical Pacific and the Indian Ocean (for the ensemble mean), as well as positive correlations 
with the North Atlantic SST; the spatial pattern of these correlations over the ocean indeed resembles 
the spatial pattern of the SST anomalies shown in Fig. 10. 

In summary, there appear to be significant temporal correlations between JJA surface temperature 
over large regions of Eurasia and SST, particularly in the North Atlantic and the tropical Pacific. 
However, the SST-forced response appears to be intertwined with and sensitive to the excitation of the 
basic internally-generated Rossby wave structures discussed previously. This aspect of the response is 
currently not well understood. 


5. Predictability and Projections 

In this section we review and provide new results on the predictability of drought and heat waves. We 
also review studies that examine longer-tenn projections, including those that examine overall trends 
in precipitation and temperature and provide an outlook for future heat waves and droughts. 
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a) A Review of Predictability 

The predictability of heat waves and associated droughts is particularly challenging in view of their 
strong link to the development of persistent anticyclones, blocking, and stationary Rossby waves (see 
previous section on mechanisms). Most weather and climate models do not adequately represent 
blocking events (e.g., Scaife et al. 2010); they underestimate the occurrence of blocking as well as its 
intensity and duration. In addition, the basic predictability of blocking is likely rather short (perhaps 
a few weeks), since blocking ridges are believed to be maintained by interactions with smaller scale 
weather systems (e.g., Scaife et al. 2010). The aforementioned results tying the development of major 
Eurasian heat waves to stationary Rossby waves also indicates relatively short predictability time 
scales, since the main forcing of such waves appears to be sub-monthly weather transients (e.g., 
Schubert et. al. 2011). We note that the link between Rossby waves and the development of blocking 
events is still unclear (e.g., Nascimento and Ambrizzi 2002; Woollings et al. 2008). 

Soil moisture anomalies and associated land-atmosphere feedbacks do provide some hope for skillful 
predictions out to perhaps 2 months (i.e., beyond weather time scales), though the levels of attainable 
skill, particularly given the observational networks available for soil moisture initialization, are 
modest at best (Koster et al. 2012; see also Volodin 2011). As discussed above, it is still an open 
question whether links to SST variability are sufficiently robust to provide useful forecast skill at 
seasonal to interannual time scales, though there is some evidence that Arctic sea ice changes could 
provide some predictable signals. Modes of variability such as the NAO, the NAM, East 
Atlantic/Western Russia (EA/WR) and the Scandinavian modes, while primarily associated with cold 
season variability, can also play a role by preconditioning the soil moisture for the subsequent 
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summer. While the NAO appears to have limited predictability on monthly and longer time scales 
(e.g., Johansson 2007), being largely driven by internal atmospheric dynamics, some evidence 
suggests that predictability may be provided through the coupling of the NAM with the stratosphere 
(e.g., Kornich 2010). Basic understanding of the mechanisms and predictability of the EA/WR and 
Scandinavian modes is not well understood. Monsoonal flows (e.g., Trenberth and Fasullo 2012) 
provide another potential source of predictability, though the current capability of simulating Asian 
monsoon variability is quite limited. 

The 2010 heat wave provides an important example of our current ability to predict a particularly 
extreme event. Matsueda (2011), using medium range ensemble forecasts, showed some success in 
predicting aspects of the blocking and extreme surface temperatures associated with the event out to a 
lead time of 9 days, though the later stages of the blocking in early August were less well predicted, 
with most models predicting a too early decay of the blocking. Ghelli et al. (2011) found signs of the 
developing heat wave about 3 weeks in advance in predictions with the ECMWF’s suite of models, 
though the full amplitude of the event was not predicted until about 1 week in advance. These results 
are consistent with the study by Dole et al. (2011), which found no change in the probability of 
prolonged blocking events over western Russia during July 2010 for forecasts initialized in early June 
of that year, compared with hindcasts initialized in early June of other years (1981-2008). 

b) A Case Study 

Flere, using a more idealized approach, we present new results concerning the basic predictability of 
extreme heat waves and associated drought events. We examine the predictability of one of the most 
extreme events to occur over European Russia in our multi-decadal GEOS-5 AMIP simulations: the 
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extreme heat event simulated by one of our ensemble members in the summer of 2001 (see discussion 
of Fig. 13). We remind the reader that the fact that this event happened to occur in 2001 in the 
simulation (rather than 2010 as in nature) appears to be purely by chance, since there is no consistency 
among the various ensemble members as to the timing of such events. Here we chose the event that 
occurred in 2001 in ensemble member 6 because it was one of the most extreme simulated events to 
occur over basically the same region as the observed event of 2010. Figure 15 shows the evolution of 
the surface air temperature, upper tropospheric meridional wind, and soil moisture from May through 
August of that year. The surface air temperature anomalies during May show a wave structure across 
northern Eurasia, the same structure that characterizes the monthly variability of observations (e.g., 
Fig. 3). At this time the largest temperature anomalies occur over eastern Siberia, with negative 
anomalies to the west and positive anomalies over Eastern Europe and European Russia. The same 
basic structure continues into June, showing some propagation to the east and intensification of the 
warm anomalies over European Russia, especially just north of the Black Sea, where it achieves its 
maximum amplitude of more than 5°C. By July the wave structure is more diffuse, but the warm 
anomalies over European Russia continue through July and well into August. The upper level wind 
shows that the Ts anomalies are associated with Rossby wave-like structures that develop in May, 
peak in June, and dissipate thereafter. The soil moisture anomalies show the same basic wave 
structure, though somewhat phase shifted to the east of the temperature anomalies. The negative soil 
moisture anomalies over European Russia are already evident in May (just north of the Black and 
Caspian Seas), intensifying in June and continuing through July into August. As the Ts anomalies 
move to the east, they appear to become phase locked with the soil moisture anomalies beginning with 
July and extending into August. 
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Our interpretation of the above results is that the heat wave was initiated by the development of a 
Rossby wave (May and June). This wave generated Ts anomalies that eventually became phase 
locked with existing dry soil moisture anomalies over European Russia, which acted to intensify and 
persist the Ts anomalies beyond the lifespan of the Rossby wave. We note that central ER 
(specifically the region 30°E-48°E, 50°N-60°N) in that ensemble member experienced soil moisture 
deficits for almost a decade from the mid 1990s to early 2003, deficits that appear to be part of a 
general drying and warming trend (evident also in the ensemble mean) that begins in the mid-20 th 
century and becomes especially pronounced after the mid-1990s. This suggests that the SST/GHG 
forcing may have set the stage for the development of the extremely warm summer over European 
Russia. We note that the evolution described here is quite similar to that found in Lyon and Dole 
(1995) for the 1980 and 1988 U.S. drought cases, where anomalous wavetrains associated with early 
stages of heat wave/droughts became very weak by early July, with reductions in evapotranspiration 
over the drought regions intensifying and prolonging the excessive heat into later summer. It remains 
to be seen whether GHGs may have set the stage for the development in these observed cases or 
whether naturally occurring drought would lead to the same outcome, 

We investigate the predictability of this event by performing a supplemental set of 20 simulations, 
each initialized on Oz, 15 May 2001 and run through August. Each simulation differs from the control 
(i.e., the ensemble member that produced the extreme event in Figure 15) only in the initialization of 
the atmosphere; to produce the atmospheric initial conditions, a small perturbation was added to the 
control atmosphere’s state on May 15. The results are presented in Figure 16. The left column shows 
the Ts anomalies from the control simulation, the middle column shows those for the ensemble mean 
of the perturbation experiments, and the right column shows the ensemble mean of the soil moisture 
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anomalies. In the ensemble mean, the wave structure in Ts is largely gone by June, indicating that in 
the span of a few weeks, the Rossby wave producing it has already lost all predictability. What 
remains in June is a general warm anomaly along the [40°N,50°N] latitudinal belt of continental 
Eurasia. This warming lasts into July and August over southern European Russia (the core of the 
original Ts anomaly) and is collocated with the dry soil wetness anomalies, suggesting that the land 
anomalies act to maintain the Ts anomalies several months beyond the predictability limit of the 
Rossby wave. 

In summary, for at least some extreme heat wave and drought events, predictability associated with 
stationary Rossby waves, which are largely forced by submonthly transients, appears to be limited to 
perhaps 2-3 weeks. Nevertheless, there appears to be some longer-term predictability tied to the 
persistence of soil moisture anomalies. Ties to SST variations could provide some predictability on 
seasonal and longer time scales, though SST impacts appear to be intertwined with the underlying 
internally forced and shorter-time scale Rossby wave structures, and this connection is currently 
poorly understood. 

c. A Review of projections for the future 

We now address the question of how heat waves and droughts might manifest themselves in a future, 
warmer world. 

Galos et al. (2007) reviewed drought occurrence in Hungary, noting that annual mean temperatures 
became warmer in the second half of the 20 th century, accompanied by a significant increase in 
drought frequency. In particular, summers for the period 1990-2004 were warmer than those of the 
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previous 30 years. The period 1983-1994 was an extraordinarily dry period, with severe droughts in 
the Carpathian Basin of Hungary. They found from an analysis of REMO 21 st Century simulations (a 
limited area model forced by lateral boundary conditions from three different ECHAM5/MPI-OM 
GCM runs - IPCC scenarios Bl, A IB, and A2) that the probability of dry summers will not increase in 
the first half of the 2 1 th Century, but the intensity of dry events will increase due to the higher 
temperatures. They also found, however, that during the second half of the 21 st century both the 
number and intensity of dry events will increase significantly. 

Meehl and Tebaldi (2004) examined simulations of 20 th and 21st Century climate produced by the 
global Parallel Climate Model (PCM), which used a “business as usual” emission scenario for the 21 st 
Century. They found the circulation patterns associated with heat waves in North America and 
Europe to be intensified in the 21 st Century, implying that future heat waves will be more intense, 
more frequent, and longer lasting in the second half of that century. 

Barriopedro et al. (2011) show that the 2003 and 2010 summer heat waves likely produced the 
warmest seasonal temperatures seen in 500 years over about 50% of Europe. They conclude, based on 
regional climate model simulations driven by different GCMs forced by A1B emission scenarios, that 
the probability of a summer mega-heat wave over Europe will increase by a factor of 5 to 10 in the 
next 40 years, though the probability of an event with the magnitude seen in 2010 will remain 
relatively low until the second half of the 21 st century. 

The Special Report of IPCC on extremes (IPCC, 2012) gives drought projections low confidence due 
insufficient agreement among the individual projections resulting from both model differences and 
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dependencies on the definition of drought (e.g., soil moisture versus precipitation-based indices). On the 
other hand, they conclude that is very likely that the length, frequency, and/or intensity of heat waves 
(defined with respect to present regional climate) will increase over most land areas. In particular over the 
high latitudes of the Northern Hemisphere, a 1 -in-20 year annual hottest day is likely to become a l-in-5 year 
annual extreme by the end of the 21st century, under the SRES A2 and A IB emission scenarios. 


6. Summary and Concluding Remarks 

Drought and heat waves often go hand in hand. While this can of course simply be because drier soils 
produce less evaporative cooling of the surface, in northern Eurasia persistent anticyclones appear to play a 
key role, acting to both warm and dry the atmosphere and land surface over many important agricultural 
regions, from European Russia to Kazakhstan and beyond. The importance of anticyclones in the 
development of droughts was kn own as far back as the early 20 th Century (e.g., Buchinsky 1976; also the 
review of earlier literature in Obukhov et al. 1984). Different air masses are linked to the development of 
anticyclones, especially the intrusion of Arctic air masses that occasionally combine with subtropical air 
(e.g., associated with the Azores high in eastern Europe and western Russia); a basic understanding for how 
these air masses produce especially severe droughts across Eurasia was already established by that time (e.g., 
see summary by Kleshenko et al. 2005). Perhaps less well understood, though mentioned in early historical 
documents, was the tendency for especially severe droughts and heat waves to be juxtaposed with wet and 
cool conditions in regions thousands of miles to the east or west (see Section 1). Observational studies also 
established that while droughts across northern Eurasia rarely last for more than 2 months (Cherenkova, 
2008), there is considerable evidence for longer (even multi-year) droughts to occur in the more southern 
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marginal semi-arid and arid regions of northern Eurasia (e.g. in Kazakhstan during the 1930s: cf., Kahan 
1989; Almaty 2006, or over the Great Steppe of Central and East Asia: cf., Gumilev 1960). 


Here we provide an updated picture of the role of anticyclones in northern Eurasian summer climate 
through an analysis of the last three decades of monthly surface temperature and precipitation 
variability and covariability, using the latest generation of reanalyses and gridded station observations. 
We also examine longer-term changes (including the recent decadal changes) in surface temperature 
and precipitation over Eurasia and the interannual variability of these quantities over the last century 
or so, using model simulations (especially those with the GEOS-5 AGCM) to better understand the 
nature of the variability. 

Among the key new results of this study is the quantification of the major summer patterns of monthly 
surface temperature and precipitation variability across northern Eurasia and the link between these 
patterns and stationary Rossby waves. The characteristic east- west wave structure of the leading 
patterns of surface meteorological variables are a reflection of these waves which, when amplified and 
stationary, appear to have led to some of the most extreme heat waves and droughts in Eurasia (e.g., 
the 2003 European and 2010 Russian heat waves), with anomalies of opposite sign occurring to the 
east and/or west depending on the phase and location of the wave. These waves appear to be initially 
forced upstream of Eurasia (e.g., within the North Atlantic jet exit region [Schubert et al. 2011]); the 
wave energy propagates over northern Eurasia, north of the mean jet and/or further to the south where 
it remains confined to the mean jet. The structure of these waves and their time scales (weeks to a 
few months) are consistent with past observations of the structure and time scales of heat waves and 
droughts across northern Eurasia. 
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The GEOS-5 AGCM simulations forced with observed SSTs and GHGs show heat waves that appear 
to be linked to Rossby waves occurring over Eurasia, including some rare, very extreme events during 
the last few decades. A case study of one of the most extreme heat waves to occur in the model 
(during the “summer of 2001” of one ensemble member) shows that the associated Rossby-like wave 
pattern in the surface temperature anomalies is for the most part unpredictable beyond about 1 month. 
Some aspects of the heat wave are, however, predictable for several months: these are the surface 
temperature anomalies at the center of the heat wave associated with soil moisture anomalies that 
persist through the summer. An inspection of the precursors to the heat wave show existing dry soil 
moisture anomalies (especially pronounced in that ensemble member) that are part of a long term 
drying and warming trend simulated in the model, a trend that is consistent with observations. More 
generally, the impact of land-atmosphere feedbacks was quantified with model simulations in which 
the soil moisture feedbacks were disabled. These runs show that temperature variability is especially 
strongly tied to soil moisture variability in the southern parts of our study area extending from 
southern Europe eastward across the Caucasus, Kazakhstan, Mongolia, and northern China. 

Our investigation of the warming that has been observed over northern Eurasia in the last three 
decades shows that it is part of a large-scale pattern of wanning with local maxima over European 
Russia and over Mongolia/eastern Siberia. Precipitation changes consist of deficits across Eurasia 
covering parts of northeastern Europe, European Russia, Kazakhstan, southeastern Siberia, Mongolia, 
and northern China. Precipitation increases occur across Siberia north of about 60°N. Remarkably, 
the ensemble mean of the AGCM simulations forced with observed SST and GHGs to a large extent 
reproduces the observed surface temperature and precipitation trend patterns of the last three decades, 
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though with smaller amplitude. This suggests that some of the basic features of the observed trends 
over Eurasia are associated with an SST trend that consists of a PDO-like colder Pacific and an AMO- 
like warmer Atlantic. Various model simulations (Schubert et al. 2009) carried out with idealized 
versions of these basic SST patterns indicate a global scale response to the PDO-like and AMO-like 
patterns, a response that is intensified by a global warming SST trend pattern. The dynamical 
response of the models to the SST forcing consists of a zonally symmetric positive upper tropospheric 
height anomaly in the middle latitudes of both hemispheres that appears to provide the large-scale 
atmospheric teleconnections linking the various regions of the world. We speculate that such a 
response was responsible for the synchronicity of droughts in such disparate regions as the Eurasian 
grainbelt (spanning Russia, Ukraine, and Kazakhstan) and the US Great Plains during, for example, 
the 1930s, as well as the drought and extreme heat in the same regions during the summer of 2012. It 
is also suggested that the longer time scales of dry conditions in the more southern regions of northern 
Eurasia may be induced by global SST anomalies. 

A survey of the literature indicates a general consensus that the future holds an enhanced probability 
of heat waves across northern Eurasia especially by the second half of the 21 st Century, while there is 
less certainty regarding future drought, reflecting the greater uncertainty in precipitation and soil 
moisture projections compared with temperature. It is also clear that there are still gaps in our 
understanding of the physical mechanisms that control the intensity, duration and frequency of heat 
waves and droughts. Perhaps most important are the uncertainties that remain in our understanding 
of the interactions between the short-term atmospheric variability associated with extremes, and the 
longer-term variability and trends associated with soil moisture feedbacks, SST anomalies, and an 
overall warming world. 
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Appendix A: Observational Datasets and Model Simulations 
Al . Observations and Reanalyses 

Our analysis is based in part on MERRA (Rienecker et al. 2011). MERRA is an atmospheric reanalysis that 
was produced with the Goddard Earth Observing System Data Assimilation System Version 5 (GEOS-5) 
documented in Rienecker et al. (2008), consisting of the GEOS-5 atmospheric model and the Grid-point 
Statistical Interpolation (GSI) analysis system, the latter being a system jointly developed by the GMAO and 
NOAA’s National Centers for Environmental Prediction. The GEOS-5 assimilation system includes an 
incremental analysis update (IAU) procedure (Bloom et al. 1996) that slowly adjusts the model states toward 
the observed state. This has the benefit of minimizing any unrealistic spin down (or spin-up) of the water 
cycle. MERRA was run at a resolution of Vi latitude x 2/3° longitude with 72-levels extending to O.OlhPa. 
More information about MERRA can be found at: http://gmao.gsfc.nasa.gov/research/merra/ . The MERRA 
data used in this study (Ts, T2m, Precip, V250) were all taken at the full resolution of 'A° latitude x 2/3° 
longitude covering the period 1979-2010. Limited comparisons are made with ECMWF’s ERA-interim 
reanalysis (Dee et al. 201 la, b). 

We also make use of various station observations. These are the NOAA NCEP CPC GHCN_CAMS 
gridded two-meter temperature (Fan et al. 2008) - a 0.5°X0.5° degree latitude/longitude resolution 
dataset covering the period January 1948-January 2013. We also make use of the CRUTEM4 two 
meter temperature station data gridded to 5° X 5° latitude/longitude for the period 1850-2012 (Jones et 
al. 2012), and the NOAA Merged Land-Ocean Surface Temperature Analysis (MLOST, Voss et al. 

2012), V3.5.2 also at 5° X 5° latitude/longitude for the period 1880-present. For the precipitation 
data, we use NOAA’s precipitation reconstruction over land (PRECL) on a 1° latitude/longitude grid 
for the period 1948-2013 (Chen et al. 2002). The other precipitation data used in the study is version 
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2.2 of the combined Global Precipitation Climatology Project (GPCP) data available on a 2.5° X 2.5° 
grid from 1979 to June 2011 (Adler et al. 2003). 

A2. The GEOS-5 Model and Simulations 

We take advantage of an ensemble of 12 AMIP-style simulations carried out with the NASA Goddard 
Earth Observing System (GEOS-5) Atmospheric General Circulation Model or AGCM (Rienecker et 
al., 2008) forced with observed SST for the period 1870 -2012. The runs were started from different 
atmospheric and land initial conditions. Ten of the twelve ensemble members were run with 
interactive aerosols, while the other two used a prescribed aerosol climatology. We have found no 
discemable difference in the basic climatology and time dependence due to the treatment of the 
aerosols, so for the purposes of this study our ensemble means are based on all twelve runs. We also 
present some results on the impact of soil moisture feedback (Section 3a). Those results are based on 
two thirty-three year simulations for 1980-2012 forced with observed SST 1 *. The first was run with 
interactive land, while the second was run with specified climatological soil moisture computed as an 
average of a previously run multi-decadal simulation. Details of the model are described next. 

The GEOS-5 AGCM employs the finite-volume dynamics of Lin (2004). This dynamical core is 
integrated with various physics packages (Bacmeister et al., 2006) under the Earth System 
Modeling Framework (Collins et al., 2005) including the Catchment Land Surface Model (Koster 
et al., 2000), and a modified form of the Relaxed Arakawa-Schubert convection scheme described 
by Moorthi and Suarez (1992). For the experiments described here we used version 2.4 of the 


11 In practice these two runs were reinitialized November 1 of each year from a previous long model simulation forced with 
observed SST. 
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AGCM. The model was run with 72 hybrid-sigma vertical levels extending to O.OlhPa, and 1° 
(about 100km) horizontal resolution on a latitude/longitude grid. 

The CO 2 consists of the time varying annual global mean values provided by IPCC/CMIP5. The 
other greenhouse gases (GHGs: CH 4 , N 2 0, CFC-1 1, CFC-12, and FICFC-22), stratospheric water 
vapor (H 2 0), and ozone (O 3 ) are relaxed to time varying zonal averages with a two-day e-folding 
time. The zonal averages of the GHGs are taken from simulations of 1950-2010 with the GEOS 
chemistry climate model (CCM, Pawson et ah, 2008), and are calibrated (bias corrected) to the 
tropospheric concentrations specified by CMIP5 (Meinshausen et ah, 2011). Stratospheric H 2 0 is 
also taken from the CCM. In both cases, GHGs and H 2 0, five-year running averages are first 
computed to reduce the influence of interannual variability in the CCM fields. Ozone is specified 
from AC&C/SPARC monthly averages (ftp-esg.ucllnl.org) from 1870-2005, and is converted to 
zonal means before interpolation onto GEOS-5 layers. For all seven gases, the relaxation fields 
have realistic latitudinal, vertical, and seasonal variations imposed on their specified trends. Two- 
day e-folding times allow the species contours to sufficiently follow planetary-scale potential 
vorticity defonnations in the stratosphere. 

Aerosols are computed using the Goddard Chemistry, Aerosol, Radiation, and Transport model 
(GOCART, Chin et al. 2002; Colarco et al. 2009) in GEOS-5. The GOCART module is run 
online within the GEOS-5 AGCM; that is, the aerosols and other tracers are radiatively interactive 
and transported consistently with the underlying hydrodynamics and physical parameterizations 
(e.g., moist convection and turbulent mixing) of the model. GOCART treats the sources, sinks, 
and chemistry of dust, sulfate, sea salt, and black and organic carbon aerosols. Aerosol species are 
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assumed to be external mixtures. Total mass of sulfate and hydrophobic and hydrophilic modes of 
carbonaceous aerosols are tracked, while for dust and sea salt the particle size distribution is 
explicitly resolved across 5 non-interacting size bins for each. 

Both dust and sea salt formulations have wind-speed dependent emission functions, while sulfate 
and carbonaceous species have emissions principally from fossil fuel combustion, biomass 
burning, and biofuel consumption, with additional biogenic sources of organic carbon. Sulfate has 
additional chemical production from oxidation of SO 2 and DMS, and we include a database of 
volcanic SO 2 emissions and injection heights. For all aerosol species, optical properties are 
primarily from the commonly used OPAC data set (Hess et al. 1998). This framework also 
includes the representation of CO tracers, which have emissions from fossil fuel, biofuel, and 
biomass burning. The online CO processes in GEOS-5 derive from Bian et al. (2007), and include 
indirect production of CO from oxidation of natural and anthropogenic non-methane 
hydrocarbons, chemical production from methane (CH4) oxidation, and losses through reaction 
with OH. 

Appendix B: Northern Eurasian Droughts and Heat Waves since 1875. 

Here we briefly review some of the key metrics that have been used to characterize drought and heat 
waves over Northern Eurasia. We also include some further information (in addition to that already 
provided in the text) on past droughts and heat waves over this region. Table B1 is a compilation of 
the droughts and heat waves that have occurred since 1875, based on various scientific publications as 
well as the popular literature (references are noted in the table). The table also includes information 
on the regions affected, and other auxiliary information (comments) of potential use to those 
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interested in investigating these events further. Years in bold indicate major droughts or heat waves, 
though it must be kept in mind that these are very subjective assessments of the relative severity of the 
various events, since they are based on differing metrics that emphasize varying aspects of 
meteorological, agricultural and hydrological droughts (and heat waves), over different time periods. 
As such, Table B1 should be considered as a convenient starting point for further investigation of the 
various droughts and heat waves that have occurred over Northern Eurasia in the last 135 years or so, 
rather than an objective comparative assessment of all droughts and heat waves that have occurred 
over this very large region. In fact we view this table (in the spirit of the GDIS effort mention earlier) 
as the starting point for a continually evolving catalog of historical droughts and heat waves that have 
occurred world- wide. 


Turning now to some of the popular metrics of drought, Ped (1975) introduced the index of aridity S; 
defined as 


AT_AP_ 


(Bl) 


o T o p 

where AT and A R are the deviations (from a long term mean) of monthly mean air temperatures and 
precipitation, and o T and o p are their standard deviations. This index has been used frequently in the 
CIS for identification of atmospheric drought in terms of three classes: light (0< S,<2.0), average (2.0< 
Si<3.0), and strong (Si>3.0). 


Another popular index is the hydrothermal coefficient (HTC) developed by Selyaninov (1928). 


HTC = 



(B2) 
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where the numerator is the total rainfall over the growing season (in mm) and the sum in the 
denominator is the accumulated mean daily surface air temperature above 10°C for the same time 
period. The threshold for drought is typically HTC values less than or equal to 0.8, with severe 
droughts having HTC values of 0.4 or less. 

Meshcherskaya and Blazhevich (1997) developed a combined drought and excessive moisture index 
(DM) that takes into account the areal extent of the precipitation and temperature anomalies. Drought 
(excessive wet) conditions are defined according to whether the precipitation falls below (exceeds) 

80% (120%) of the long term mean, and the temperature anomalies are above (below) 1°C (-1°C). 

They produced a catalog of drought occurrence over the main grain producing regions of the FSU for 
May-July 1891-1995 for both the European and Asian parts (see their Table 4). They found that the 
most severe droughts (in order of decreasing severity) in the European region occurred during 1936, 
1975, 1979, and 1891, while for the Asian part the most severe droughts occurred during 1955, 1965, 
1951, and 1931. 

Dai (2011) compared different forms of the PDSI, finding generally little difference between four 
different formulations. They generally compare well with monthly soil moisture observations 
(Robock et al 2000), annual streamflow and monthly GRACE satellite-observed water storage 
changes. For example, correlations of up to 0.77 were found in parts of the FSU for soil moisture in 
the top lm even over high latitude cold regions (east of the Urals). 

Another more recently developed drought index that includes the effects of temperature on drought 
variability is the standardized precipitation evapotranspiration index (SPEI, Vicente-Serrano et al. 2010). 
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The SPEI, similar to the sc-PDSI (Wells et al. 2004), can capture increases in drought severity associated 
with higher water demand as a result of evapotranspiration, under global warming conditions. The SPEI was 
used for example by Potop and Mozny (2011) to study the evolution of drought in the Czech Republic. They 
found that increasing temperatures played a role in the intensification of the droughts during the 1980s and 
1990s. 

Rocheva (2012) proposed a 500mb height index as an indicator of drought over the main grain- 
producing regions of Russia using NCEP/NCAR reanalysis data. They found that droughts occurred 
during 1972, 1975, 1979, 1981, 1995, and 1998 (they note that their findings are consistent with the 
droughts identified by Strashnayaa and Bogomolova, 2005). 

According to Golubev and Dronin (2004), droughts in Russia during the last hundred years tended to 
occur over three main geographical areas consisting of Central, Southern and Eastern Russia (based 
on the work of Anon, 1933). The central type of drought covers the Volga basin, the Northern 
Caucasus and the central chernozem region and some oblasts of the central region, affecting the major 
agricultural regions of Russia, and the forest zone of European Russia (associated with numerous 
forest fires in the central and northern regions). The southern type of drought is limited to the Volga 
basin and Urals region, and while it covers less area, its intensity has generally been more severe and 
has often destroyed the entire crop production of the affected region. They note that the eastern type of 
drought affects the steppe and forest-steppe of Siberia and this usually occurs when the southern part 
of European Russia is characterized by good weather. This again highlights the juxtaposition of 
drought and wet conditions as a characteristic feature of climate variability over Eurasia - in this case 
the contrast is between European and Asian Russia. Golubev and Dronin (2004) summarize the past 
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occurrence of each type of drought with the central droughts occurring during 1920, 1924, 1936, 
1946, 1972, 1979, 1981 and 1984, the southern droughts occurring during 1901, 1906, 1921, 1939, 
1948, 1951, 1957, 1975, 1995, and the eastern droughts occurring during 1911, 1931, 1963, 1965, 
1991. One of the worst modem droughts over ER occurred in the summer of 1972 (Fedorov 1973; 
Buchinsky 1976). That drought was associated with an anticyclone that was centered over Moscow 
and that established in May and persisted throughout the summer. The drought appears to have 
started in eastern Ukraine and was at the time characterized as a 100-year event. 

A NOAA team of experts 12 note that western Russia has a climatological vulnerability to blocking (see 
also Tyrlis and Hoskins 2008; Woollings et al. 2008; Dole et al. 2011) and associated heat waves (e.g., 
1960, 1972, 1988). They point out, however, that a high index value for blocking days is not a 
necessary condition for high July surface temperature over western Russia— e.g., the warm summers 
of 1981, 1999, 2001, and 2002 did not experience an unusual number of blocking days. 

Almaty (2006) found that hydrological drought (low runoff) occurred in the western regions of 
Kazakhstan during 1933-1939, 1972-78 and 1996-97. The latter two periods were also low-water 
periods in northern Kazakhstan, while 1963-65, 1967-70 and 1996-2000 were low-water periods in 
central Kazakhstan. Drought comes to the lowland of southern Kazakhstan roughly every 4-5 years. 

It was in drought during 2000-2001 with the Chardarya reservoir having the lowest water storage 
since 1977 in August of 2001. 


12 http://www.esrl.noaa.gov/psd/csi/events/2010/russianheatwave/prelim.html 
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1732 Table B 1: Droughts and heat wave years in Eurasia ( 1875-2012). ENSO phase refers to summer 

1733 months. Bold indicates major droughts/ heat waves as identified by various authors. 


Year 

Months 

Regions Affected 

Source (s) 

Comments 

1875 


Ukraine, Non-Black Soil Region, North Caucasus 

Kahan 1989; Rudenko 1958 

General Dmitry Milyutin’s 
diary (see text); Strong La 
Nina 

1879 


Drought in middle Volga, central Ukraine 

Kahan 1989 

La Nina 

1885 


South and East Ukraine, middle and lower Volga, North 
Caucasus 

Kahan 1989 

Tchaikovsky letter 13 

1889 


Drought in Ukraine, Lower Volga 

Kahan 1989 

La Nina 

1890 


Drought in Central and southern regions of European Russia, 
Lower Volga, Ural; severe drought in Ukraine 

Fedorov 1973; Buchinsky 
1976; Kahan 1989; Boken et 
al. 2005 

La Nina 

1891 

May- 

August 

Drought in Central and southern regions of European Russia, 
spreading north to middle course of Kama and Vyatka, 
Blacksoil (Chernozem) Region, all of Volga, Southern Ukraine, 
North Caucasus, 80% of European Territory of USSR (ETU) 
between 45-55N; 30-50E affected by drought 

Fedorov 1973; Buchinsky 
1976; Meshcherskaya and 
Blazhevich 1997; Kahan 
1989; Polozova and 
Grigoryeva 1984 

Absolute temperature record 
for May was set (+31.8 deg. 
C), not broken until 2007 

1892 


Drought in Central and southern regions of European Russia; 
Ukraine, Central Blacksoil Region, Lower Volga, North 
Caucasus; 90% of (ETU) between 45-55N; 30-50E affected by 
drought (Polozova and Grigoryeva 1984) 

Fedorov 1973; Buchinsky 
1976; Kahan 1989; Polozova 
and Grigoryeva 1984 

La Nina 

1897 


Drought in South Ukraine, Lower Volga; 80% of (ETU) 
between 45-55N; 30-50E affected by drought 

Kahan 1989; Polozova and 
Grigoryeva 1984 


1898 


Severe drought in Ukraine 

Boken et al 2005 


1900 


Drought in Ukraine, Western Siberia 

Kahan 1989 

El Nino 

1901 


“Southern” type of drought in the Volga River Basin (Volga 
and Volga- Vyatka) and the Urals Region; eastern Ukraine, 
Ural, lower and middle Volga; 95% of (ETU) between 45-5 5N; 
30-50E affected by drought 

Golubev and Dronin, 2004; 
Kahan 1989;Polozova and 
Grigoryeva 1984 


1905 


60% of (ETU) between 45-55N; 30-50E affected by drought 

Polozova and Grigoryeva 
1984 


1906 


“Southern” type of drought in the Volga River Basin (Volga 
and Volga- Vyatka) and the Urals Region; eastern Ukraine, 
Volga; 75% of (ETU) between 45-55N; 30-50E affected by 
drought 

Golubev and Dronin, 2004; 
Kahan 1989; Polozova and 
Grigoryeva 1984 

Severe drought in Ukraine 

1911 


“Eastern” type of drought in-steppe and forest-steppe of 
Western and Eastern Siberia; European Russia, Ukraine, west 
Siberia; severe drought in Ukraine; 65% of (ETU) between 45- 
55N; 30-50E affected by drought 

Golubev and Dronin, 2004; 
Kahan 1989; Boken et al 
2005; Polozova and 
Grigoryeva 1984 

La Nina 

1914 


Volga, Ukraine; 75% of (ETU) between 45-55N; 30-50E 
affected by drought 

Kahan 1989; Polozova and 
Grigoryeva 1984 

El Nino 

1917 


65% of (ETU) between 45-55N; 30-50E affected by drought 

Polozova and Grigoryeva 
1984 


1918 


Ukraine 

Rudenko 1958 


1920 


Drought in Central Russia, forest zone of European Russia; 
Ukraine, Volga; 95% of (ETU) between 45-55N; 30-50E 
affected by drought 

Golubev and Dronin, 2004; 
Kahan 1989; Polozova and 
Grigoryeva 1984 


1921 


“Southern” type of drought in the Volga River Basin (Volga 
and Volga-Vyatka) and the Urals Region; European Russia, 
Ukraine, Volga, Western Siberia; 95% of (ETU) between 45- 
55N; 30-50E affected by drought 

Golubev and Dronin, 2004; 
Kahan 1989; Rudenko 1958; 
Polozova and Grigoryeva 
1984 

1921 Famine in Russia 
(Volga-Ural region); La Nina; 
severe drought in Ukraine 

1924 


Drought in Central Russia, forest zone of European Russia; 
Central Blacksoil Region, Lower Volga; 90% of (ETU) 
between 45-55N; 30-50E affected by drought 

Golubev and Dronin, 2004; 
Kahan 1989; Polozova and 
Grigoryeva 1984 

La Nina 

1931 

July 

Heat wave; “Eastern” type of drought in steppe and forest- 
steppe of Western and Eastern Siberia; Ural, Central Blacksoil 
Region; 75% of (ETU) between 45-55N; 30-50E affected by 
drought 

NOAA expert team 14 ; 
Golubev and Dronin, 2004; 
Kahan 1989; Meshcherskaya 
and Blazhevich 1997; 
Polozova and Grigoryeva 



13 (in a letter from Tchaikovsky (composer) to N. von Mekk, 3-10 August 1885, in Majdanovo, now Klin): “I’m writing to you at 
three o’clock in the afternoon in such darkness, you would think it was nine o’clock at night. For several days, the horizon has 
been enveloped in a smoke haze, arising, they say , from fires in the forest and peat bogs. Visibility is diminishing by the day, and 
I’m starting to fear that we might even die of suffocation.” http:lltherese-phil.liveiournal.comll71196.html 

14 http://www.esrl.noaa.gov/psd/csi/events/2010/russianheatwave/prelim.html 
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1984 


1933-39 


Drought (low-water) years in western Kazakhstan 

Almaty 2006 


1934 


70% of (ETU) between 45-55N; 30-50E affected by drought 

Polozova and Grigoryeva 
1984 


1936 


Drought in Central Russia, forest zone of European Russia; 
Ukraine, Volga, Western Siberia; 90% of (ETU) between 45- 
55N; 30-50E affected by drought 

Golubev and Dronin, 2004; 
Meshcherskaya and 
Blazhevich 1997; Kahan 
1989; Polozova and 
Grigoryeva 1984 

Dry preceding winter and 
spring - Cherenkova 2012 

1938 


Central Blacksoil Region, Eastern Ukraine; 85% of (ETU) 
between 45-55N; 30-50E affected by drought 

Kahan 1989; Polozova and 
Grigoryeva 1984 

Dry preceding winter and 
spring - Cherenkova 2012; La 
Nina 

1939 


“Southern” type of drought in the Volga River Basin (Volga 
and Volga-Vyatka) and the Urals Region; Lower Volga, Upper 
Volga, Southeast Ukraine; 65% of (ETU) between 45-55N; 30- 
50E affected by drought 

Golubev and Dronin, 2004; 
Kahan 1989; Polozova and 
Grigoryeva 1984 


1946 


Drought in Central Russia, forest zone of European Russia; 
Ukraine, Central Blacksoil Region, Volga, North Caucasus; 
100% of (ETU) between 45-55N; 30-50E affected by drought 

Golubev and Dronin, 2004; 
Kahan 1989; Polozova and 
Grigoryeva 1984 

Famine; serious crop failures 
in grain production in USSR; 
drought covered 50% of total 
agricultural land of USSR 

1948 


“Southern” type of drought in the Volga River Basin (Volga 
and Volga-Vyatka) and the Urals Region 

Golubev and Dronin, 2004; 
Kahan 1989 

serious crop failures in grain 
production in USSR 

1950 


Eastern Ukraine, Middle and Lower Volga; 85% of (ETU) 
between 45-55N; 30-50E affected by drought 

Kahan 1989; Polozova and 
Grigoryeva 1984 

La Nina 

1951 


“Southern” type of drought in the Volga River Basin (Volga 
and Volga-Vyatka) and the Urals Region; drought in Asian part 
of FSU; Ukraine, Lower Volga; 75% of (ETU) between 45- 
55N; 30-50E affected by drought 

Golubev and Dronin, 2004; 
Meshcherskaya and 
Blazhevich 1997; Kahan 
1989; Polozova and 
Grigoryeva 1984 

El Nino 

1953 


Ukraine, Volga 

Kahan 1989 


1954 


Southern Ukraine, Lower Volga; 60% of (ETU) between 45- 
55N; 30-50E affected by drought 

Kahan 1989; Polozova and 
Grigoryeva 1984 

La Nina 

1955 

July 

Heat wave over Russia; drought in Asian part of FSU; western 
Siberia, Lower Volga 

Meshcherskaya and 
Blazhevich 1997; Kahan 
1989 

La Nina 

1957 


“Southern” type of drought in the Volga River Basin (Volga 
and Volga-Vyatka) and the Urals Region; 75% of (ETU) 
between 45-55N; 30-50E affected by drought 

Golubev and Dronin, 2004; 
Polozova and Grigoryeva 
1984 

serious crop failures in grain 
production in USSR; El Nino 

1959 


60% of (ETU) between 45-55N; 30-50E affected by drought 

Polozova and Grigoryeva 
1984 

serious crop failures in grain 
production in USSR 

1960 


Heat wave in Western Russia 

NOAA expert team 

Associated with blocking 

1963 


Eastern” type of drought in steppe and forest-steppe of Western 
and Eastern Siberia; western Siberia, Kazakhstan, Ural, eastern 
Ukraine; 95% of (ETU) between 45-55N; 30-50E affected by 
drought 

Golubev and Dronin, 2004; 
Kahan 1989; Polozova and 
Grigoryeva 1984 

serious crop failures in grain 
production in USSR; El Nino 

1965 


Eastern” type of drought in steppe and forest-steppe of Western 
and Eastern Siberia; drought in Asian part of FSU; western 
Siberia, Kazakhstan, central Asia, Ural, Volga 

Golubev and Dronin, 2004; 
Meshcherskaya and 
Blazhevich 1997; Kahan 
1989 

El Nino 

1972-78 


Drought (low-water) years in western and northern Kazakhstan 

Almaty 2006 


1972 

Established 
in May and 
lasted 
through 
summer 

Drought and heat wave: Anticyclone centered over Moscow, 
covered western Russia, drought started in eastern Ukraine; 
100% of (ETU) between 45-55N; 30-50E affected by drought 

Fedorov 1973; Buchinsky 
1976; Strashnaya and 
Bogomolova, 2005; 
Polozova and Grigoryeva 
1984; Cherenkova 2012 

Blocking, one of worst 
modern droughts - “100-year” 
drought; El Nino; serious crop 
failures in grain production in 
USSR; Dry preceding winter 
and spring. 

1975 


“Southern” type of drought in the Volga River Basin (Volga 
and Volga-Vyatka) and the Urals Region 

Golubev and Dronin, 2004; 
Meshcherskaya and 
Blazhevich 1997; Strashnaya 
and Bogomolova, 2005 

serious crop failures in grain 
production in USSR; La Nina 

1979 


Drought in Central Russia, forest zone of European Russia 

Golubev and Dronin, 2004; 
Meshcherskaya and 
Blazhevich 1997; Strashnaya 
and Bogomolova, 2005 

serious crop failures in grain 
production in USSR 

1981 

July heat 
wave 

Drought and heat wave in Central Russia, forest zone of 
European Russia 

Golubev and Dronin, 2004; 
Meshcherskaya and 
Blazhevich 1997; Strashnaya 
and Bogomolova, 2005; 
NOAA expert team 

No unusual blocking; 

Worst drought in European 
Russia between 1891-1995; 
serious crop failures in grain 
production in USSR 
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1734 

1984 


Drought in Central Russia, forest zone of European Russia, 
North and west UK 

Golubev and Dronin, 2004, 
Ben Lloyd-Hughes, Benfleld 
Hazard Research Centre, 
UCL 

UK: Very dry spring and 
summer led to the imposition 
of hosepipe bans. 

1988 

July 

Russia 

NOAA expert team 

Blocking, La Nina 

1991 


“Eastern 44 type drought in steppe and forest-steppe of Western 
and Eastern Siberia 

Golubev and Dronin, 2004; 

El Nino 

1992 


Germany, Hungary, Bulgaria and much of western Russia 

Ben Lloyd-Hughes, Benfleld 
Hazard Research Centre, 
UCL 

German crop production 
reduced by 22%. Irrigation 
suspended in Bulgaria. Worst 
Russian drought in 1 0 years 

1995 


“Southern” type of drought in the Volga basin (Volga and 
Volga-Vyatka) and the Urals Region 

Golubev and Dronin, 2004; 
Strashnaya and Bogomolova, 
2005 

El Nino 

1996-97 


Drought (low-water) years in western and northern Kazakhstan 

Almaty 2006 


1998 


The Volga Region, North Caucasus and Central-Chernozem, 
warmest year of the last century for Mongolia 

Strashnaya and Bogomolova, 
2005, Batima et al. 2005 


1999 - 

2002 


Summer drought in Mongolia 

Batima et al. 2005 

1999, 2000 La Nina; 
2002 El Nino 

1999 

July 

Heat wave in Western Russia 

NOAA expert team 

No unusual blocking over 
Russia; La Nina 

2001 


July heat wave in Western Russia, drought in southern 
Kazakstan 

NOAA expert team, Almaty 
2006 

No unusual blocking over 
Russia 

2002 

July 

July heat wave in Western Russia, drought in southern 
Kazakhstan 

NOAA expert team, Almaty 
2006 

No unusual blocking over 
Russia 

2007 

May, June, 
continued to 
September 
in Japan 

Heat wave in Southeast Europe (June-August), South Asian 
countries of India, Pakistan, Bangladesh and Nepal, as well as 
Russia, Japan and China 

.Buranov, Ivan, Geroyeva, 
Anna, Kornysheva, Alyona 
(29 Mav 2007). "Moscow 
Swelters in Heat Wave", 
Kommersant.com. 

The temperature in Moscow 
reached +32.9 °C (91.2 °F) on 

27 May. The Russian capital 
had not seen such a sustained 
heat wave for 128 years. On 

28 May an absolute 

temperature record for May 
was set, breaking the record 
of +31 .8 °C (89.2 °F) that had 
been registered back in 1891 

2010 

July 

Eastern Europe, Russia heat wave 

Barriopedro et al. 2011; 
Meshcherskaya et al. 2011; 
Mokhov 2011 

Pakistan: all-time high of 53.7 
C at Mohenjo-Daro, Sindh, 
May 26; Cold preceding 
winter in ETR (European 
Territory of Russia)- 
Cherenkova 2012 

2012 

summer 

Drought and heat wave over Russia, Kazakhstan, Ukraine 

http : // theextinc ti onprotoc o 1 . w 
ordpress. com/20 1 2/07/3 1/sev 
ere-heatwave-drought-and- 
wildfires-destroy-russian- 
harvest/ 

Eurasian grain belt hit hard 
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List of Figures 

Figure 1: Left panels: The mean June-July-August (JJA) precipitation variance from 

MERRA, ERA-Interim, and NOAA-PRECL observations, computed for the period 
1979-2012. Units are (mm/day) 2 . A linear trend was removed for each calendar 
month before computing the variances. Right panels: The same as the left, except 
for the variance of the Ped (1975) drought index (see text). 

Figure 2: Left panels: The mean JJA 2-meter temperature variance from MERRA, ERA- 

Interim, and GHCN CAMS observations, computed for the period 1979-2012. Units 
are (°C) 2 . A linear trend was removed for each calendar month before computing the 
variances. Right panels: The same as the left, except for the percent of the 2-meter 
temperature variance explained by precipitation (see text). 

Figure 3: The five leading Rotated Empirical Orthogonal Functions (REOFs) of the 

combined monthly fields of 2-meter temperature and precipitation for JJA of 1979- 
2012. The fields are normalized by their respective variances (standard deviation) 
and a linear trend was removed for each calendar month before computing the 
REOFs. The results are based on MERRA. Units are arbitrary. 

Figure 4: The projection of the monthly JJA normalized 2-meter temperature and 

precipitation fields onto the leading REOFs shown in Figure 3. The fields are not 
detrended. The values are referred to in the text as the Rotated Principal 
Components (RPCs). Units are arbitrary. 


82 


Figure 5: The standard deviation of the monthly JJA 2-meter temperature for the period 

1980-2012. Units are °C. The top row is for MERRA, the second row is for GEOS- 
5 AGCM simulations with interactive land, the third row is for GEOS-5 AGCM 
simulations with disabled land-atmosphere feedback, the bottom row shows the 
difference between the second row and the third row, and indicates the effect of land- 
atmosphere. 

Figure 6: The temporal correlation (shading) between the monthly JJA 250mb meridional 
wind (v) and leading RPCs of combined 2-meter temperature and precipitation for 
the period 1979-2012. A linear trend was removed for each calendar month before 
computing the covariances. Contours are the long-term mean JJA zonal wind (u) at 
250mb (15, 20, 25 m/sec). Results are for MERRA. 

Figure 7: Examples from MERRA of the 250mb v-wind anomalies (with respect to the 

1980-2010 climatology) for selected months to highlight the predilection for Rossby 
wave structures similar to those shown in Figure 6 that are linked to monthly 
variability in surface meteorology. Units: m/sec. 

Figure 8: The time evolution of the response of the eddy v-wind at c=0.257 to an idealized 
constant vorticity source at 0°E, 50°N. The results are from a stationary wave model 
(Ting and Yu 1998) with a 3-dimensional JJA mean base state taken from MERRA 
for the period 1979-2010. The horizontal structure of the idealized forcing has a 
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sine-squared functional form, with horizontal scales of 10°longitude by 10°latitude 
and vertical profile (maximum of 5.6 x e" 10 sec in the upper troposphere) following 
Liu et al (1998). See Schubert et al. (2011) for additional results. Contours are the 
long-term mean JJA zonal wind (u) at 250mb (15, 20, 25 m/sec) from MERRA. 

Figure 9: The difference between the long term JJA means (1996-2011 minus 1980-1995) for 2- 
meter temperature (°C, left panels) and precipitation (mm/day, right panels). The top row 
is for MERRA, the second row is for ERA Interim, the third row is from GFICNCAMS 
(for 2-meter temperature) and GPCP (for precipitation), and the bottom row is for the 
ensemble mean of the AMIP runs. Note the different shading intervals for the ensemble 
means. 

Figure 10: Left panel: The linear trend in the annual mean SST from FladlSST vl (Rayner 
et al 2003) for the period 1980-2011. Right panel: The idealized SST forcing pattern 
that was used in the USCLIVAR drought working group to force various climate 
models. The pattern is composed of the three leading REOFs of the annual mean 
SST consisting of the cold phase of a Pacific decadal mode, the warm phase of an 
Atlantic Muli-decadal Oscillation (AMO) - like mode, and the trend pattern (see 
Schubert et al. 2009). 


Figure 11: The JJA responses to the idealized SST forcing pattern shown in the right panel 
of Figure 10 averaged over four different AGCMs (CCM3, GEOS-5, GFS, and 
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GFDL). The top left panel is surface temperature (units: °C), the top right panel is 
for precipitation (units: mm/day), and the bottom panel is for the 200mb height 
(units: meters). 

Figure 12: The time series of the 2-meter temperature anomalies for four different regions 
across Eurasia. Top panels: Europe (10°W-30°E, 35°N-70°N); second from top: 
European Russia (25°E-60°E, 46°N-62°N); second from bottom: South-central 
Siberia (90°E-120°E, 45°N-65°N); bottom panels: Aral Sea region (45°E-75°E, 
35°N-55°N). The left panels are from the CRUTEM4 data (1871-2010). The right 
panels are based on the NOAA MLOST data (1880-2010). Units are °C. 

Figure 13: The same as Figure 12, except for two of the 12 ensemble members of the 

GEOS-5 AGCM simulations forced with observed SST and GHG forcings for the 
period 1871-2010. 

Figure 14: The correlation between the JJA mean surface temperature averaged over 

European Russia (20°E-45°E, 52°N-65°N, indicated by the box), and the surface 
temperature everywhere for the period 1901-1980. All data have a linear trend 
removed at each grid point before computing the correlations. The top panel is for 
the CRU TS3.0 observations. The middle panel is from the 12 GEOS-5 AMIP 
simulations. Here the correlations are computed for each ensemble member 
separately and then the correlations are averaged. The bottom panel is again for the 
12 AMIP simulations, but here the correlations are based on the ensemble mean. 
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Figure 15: Results from one of the most extreme heat waves in European Russia found in 
the GEOSS-5 AGCM simulations (see text for details). The left panels show the 
evolution of the 2-meter temperature anomalies (°C) from the simulation for May 
through August of 2001 (anomalies are computed with respect to the 1980-2010 
mean). The middle panels show the evolution of the 250mb v-wind (m/sec), and the 
right panels show the evolution of the surface soil wetness (dimensionless). 

Figure 16: Results from an ensemble of 20 perturbation experiments initialized on 0Z May 
15, 2001 for the ensemble member shown in Figure 15. The left panels show the 
original evolution of the 2-meter temperature (°C) from June through August of 2001 
(a repeat of part of Fig. 15). The middle panels are the same as the left panels except 
for the ensemble mean of the perturbation runs, and the right panels are the same as 
the middle panels except for the ensemble mean surface soil wetness 
(dimensionless). 
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Figure 1: Left panels: The mean June-July-August (JJA) precipitation variance from MERRA, ERA-Interim, and NOAA-PRECL 
observations, computed for the period 1979-2012. Units are (mm/day) 2 . A linear trend was removed for each calendar 
month before computing the variances. Right panels: The same as the left, except for the variance of the Ped (1975) 
drought index (see text). 
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Var(T)*Corr A 2(P,T): Observations 
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Figure 2: Left panels: The mean JJA 2-meter temperature variance from MERRA, ERA-Interim, and GHCN CAMS 
observations, computed for the period 1979-2012. Units are (°C) 2 . A linear trend was removed for each calendar 
month before computing the variances. Right panels: The same as the left, except for the percent of the 2-meter 
temperature variance explained by precipitation (see text) 
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Figure 3: The five leading Rotated Empirical Orthogonal Functions (REOFs) of the combined monthly fields of 2-meter 
temperature and precipitation for JJA of 1979-2012. The fields are normalized by their respective variances (standard 
deviation) and a linear trend was removed for each calendar month before computing the REOFs. The results are based 
on MERRA. Units are arbitrary. 
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Figure 4: The projection of the monthly JJA normalized 2-meter temperature and precipitation fields onto the leading 
REOFs shown in Figure 3. The fields are not detrended. The values are referred to in the text as the Rotated 
Principal Components (RPCs). Units are arbitrary. 
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Figure 5: The standard deviation of the monthly JJA 2-meter temperature for the period 1980-2012. Units are 
°C. The top row is for MERRA, the second row is for GEOS-5 AGCM simulations with interactive land, the third 
row is for GEOS-5 AGCM simulations with disabled land-atmosphere feedback, the bottom row shows the 
difference between the second row and the third row, and indicates the effect of land-atmosphere. 
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Figure 6: The temporal correlation (shading) between the monthly JJA 250mb meridional wind (v) and leading RPCs 
of combined 2-meter temperature and precipitation for the period 1979-2012. A linear trend was removed for 
each calendar month before computing the covariances. Contours are the long-term mean JJA zonal wind (u) at 
250mb (15, 20, 25 m/sec). Results are for MERRA. 
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Figure 7: Examples from MERRA of the 250mb v-wind anomalies (with respect to the 1980-2010 climatology) for 
selected months to highlight the predilection for Rossby wave structures similar to those shown in Figure 6 that are 
linked to monthly variability in surface meteorology. Units: m/sec. 
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Figure 8: The time evolution of the response of the eddy v-wind at o=0.257 to an idealized constant vorticity source at 
0°E, 50°N. The results are from a stationary wave model (Ting and Yu 1998) with a 3-dimensional JJA mean base state 
taken from MERRA for the period 1979-2010. The horizontal structure of the idealized forcing has a sine-squared 
functional form, with horizontal scales of 10°longitude by 10°latitude and vertical profile (maximum of 5.6 x e 10 sec 1 in 
the upper troposphere) following Liu et al (1998). See Schubert et al. (2011) for additional results. Contours are the 
long-term mean JJA zonal wind (u) at 250mb (15, 20, 25 m/sec) from MERRA 
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Figure 9: The difference between the long term JJA means (1996-2011 minus 1980-1995) for 2-meter temperature 
(°C, left panels) and precipitation (mm/day, right panels). The top row is for MERRA, the second row is for ERA 
Interim, the third row is from GHCN_CAMS (for 2-meter temperature) and GPCP (for precipitation), and the bottom 
row is for the ensemble mean of the AMIP runs. Note the different shading intervals for the ensemble means. 
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Figure 10: Left panel: The linear trend in the annual mean SST from HadlSST vl (Rayneret al 2003) for the period 
1980-2011. Right panel: The idealized SST forcing pattern that was used in the USCLIVAR drought working group to 
force various climate models. The pattern is composed of the three leading REOFs of the annual mean SST 
consisting of the cold phase of a Pacific decadal mode, the warm phase of an Atlantic Muli-decadal Oscillation 
(AMO) - like mode, and the trend pattern (see Schubert et al. 2009). 
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Figure 11 : The JJA responses to the idealized SST forcing pattern shown in the right panel of Figure 10 averaged over 
four different AGCMs (CCM3, GEOS-5, GFS, and GFDL). The top left panel is surface temperature (units: °C), the top 
right panel is for precipitation (units: mm/day), and the bottom panel is for the 200mb height (units: meters). 
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Figure 12: The time series of the 2-meter temperature anomalies for four different regions across Eurasia. Top 
panels: Europe (10°W-30°E, 35°N-70°N); second from top: European Russia (25°E-60°E, 46°N-62°N); second from 
bottom: South-central Siberia (90°E-120°E, 45°N-65°N); bottom panels: Aral Sea region (45°E-75°E, 35°N-55°N). The 
left panels are from the CRUTEM4 data (1871-2010). The right panels are based on the NOAA MLOST data 
(1880-2010). Units are °C. 
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Figure 13: The same as Figure 12, except for two of the 12 ensemble members of the GEOS-5 AGCM simulations 
forced with observed SST and GFIG forcings for the period 1871-2010. 
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Figure 14: The correlation between the JJA mean surface temperature averaged over European Russia (20°E-45°E, 
52°N-65°N, indicated by the box), and the surface temperature everywhere for the period 1901-1980. All data 
have a linear trend removed at each grid point before computing the correlations. The top panel is for the CRU 
TS3.0 observations. The middle panel is from the 12 GEOS-5 AMIP simulations. Here the correlations are 
computed for each ensemble member separately and then the correlations are averaged. The bottom panel is 
again for the 12 AMIP simulations, but here the correlations are based on the ensemble mean. 
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Figure 15: Results from one of the most extreme heat waves in European Russia found in the GEOSS-5 AGCM 
simulations (see text for details). The left panels show the evolution of the 2-meter temperature anomalies (°C) from 
the simulation for May through August of 2001 (anomalies are computed with respect to the 1980-2010 mean). The 
middle panels show the evolution of the 250mb v-wind (m/sec), and the right panels show the evolution of the surface 
soil wetness (dimensionless). 
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Figure 16: Results from an ensemble of 20 perturbation experiments initialized on 0Z May 15, 2001 for the ensemble 
member shown in Figure 15. The left panels show the original evolution of the 2-meter temperature (°C) from June 
through August of 2001 (a repeat of part of Fig. 15). The middle panels are the same as the left panels except for the 
ensemble mean of the perturbation runs, and the right panels are the same as the middle panels except for the 
ensemble mean surface soil wetness (dimensionless). 
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